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1.  OBJECTIVE 


1.1  Introduction 

Metallic  elements  volatilized  during  meteor  entry  into  the  Earth's  upper 
atmosphere  play  a  significant  role  in  the  structure  of  the  D  and  E  regions  of. 
the  ionosphere,  and,  at  least  in  the  case  of  sodium,  the  visible  day  and 
nightglow  emissions  from  the  mesosphere  and  lover  thermosphere.  Recently, 
it  was  suggested  that  sodium  and  other  meteor  metals  may  be  important  in 
stratospheric  chemistry  by  affecting  ozone  reduction  by  the  catalytic  chlorine 
cycle.6-8 

The  influx  of  meteor  metals  into  the  upper  atmosphere  has  been 

o  6—1 

estimated  to  be  3.5  x  10  kg  yr  ,  with  a  sodium  abundance  of  2Z  leading  to  a 
calculated  sodium  flux  of  1.2  x  lO'*  atoms  cm-2  s'1.  Other  estimates  of  sodium 
flux  run  as  high  as  2  x  10  cm  ‘  8  .  The  flux  of  other  metallic  species 

such  as  Mg,  Ca,  Al,  Si  and  Fe  will  be  as  much  as  10  times  higher  and 
speculative  concerns  about  their  Influence  on  upper  atmospheric  homogeneous 
and  heterogeneous  chemistry  have  been  published. 1-2» 7-8» 10 

Unfortunately,  all  attempts  to  model  the  role  of  volatilized  meteor 
metals  (particularly  sodium)  in  the  mesosphere  and  stratosphere 3” l*» 7> 1?“ 15 
have  suffered  from  an  almost  total  lack  of  measured  rate  constants.  All  such 
models  start  with  the  oxidation  of  sodium  or  other  metallic  species  in 
reaction  with  atmospheric  0,  O2,  or  03.  Until  1982,  the  only  measured 
chemical  rate  constants  available  for  any  meteor  metal  oxidation  reactions 

were  the  three  body  recombination  reactions  of  alkali  atoms  (Na,  K)  with 

n  16-18 
°2 

Na  or  K  +  0.  +  M  +  Na0„  or  K0„  +  M  .  (1) 

2  2  2 
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Flame  studies  by  Hynes  et  al 1  and  direct  measurements  at  lower 
temperatures  by  Husain  and  Plane16  showed  that  this  reaction  rate  Is  1000 
times  faster  than  previously  believed. 19-20  Recent  direct  flow  tube 
measurements18  in  our  laboratory  have  extended  these  results  from  700  K  down 
to  300  K,  and  show  that  the  reaction  rates  vary  inversely  with  temperature, 
and  confirm  the  larger  values  for  the  rate  constants. 

As  an  Indication  of  the  effect  that  the  availability  of  chemical  rate 
data  has  on  aeronomic  modeling  studies,  consider  the  theoretical  calculation 
of  several  important  sodium  oxidation  reaction  rate  constants: 

Na  +  03  -*>  NaO  +  02  (2) 

NaO  +  0  -*■  Na  +  02  (3) 

and 

NaO  +  0j  ♦  Products  (4) 

21 

performed  at  Aerodyne  Research,  Inc.  The  adoption  of  realistic  values  for 
just  these  three  rate  constants  has  led  to  a  major  revision  in  the 
understanding  of  the  structure  of  the  mesospheric  sodium  layer,  the 
magnitude  of  the  sodium  D  line  nlghtglow, 5» 21  and  the  understanding  of 
long-lived  luminous  meteor  trails.8*21 

From  the  recent  modeling  work7-8* 10» 13 » 21-22  and  comparison  of 
these  models  with  atmospheric  measurements,12*11*-15  neutral  sodium  is 
believed  to  be  transformed  via  a  series  of  chemical  reactions  involving  NaO, 
Na02,  and  NaOH.  A  survey  of  the  relevant  literature 6-8» i0* i3  provides  a 
fairly  complete  list  of  possible  neutral  sodium  reactions  of  importance. 
These  are  listed  in  Table  1.1,  in  four  groups,  which  to  some  extent,  reflect 
the  sequence  of  events  from  the  introduction  of  neutral  sodium  at  110  km 
altitude  to  its  disappearance  near  70  km.  Also  listed  are  reaction 
anthalpies,  with  associated  uncertainties. 
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Table  1.1  -  Atmospheric  Reactions  of  Neutral  Sodium  Species 


-AH°r(kJ  mol-1) 


The  main  removal  mechanisms  of  sodium  are  by  reaction  with  O2  or  03,  the 
latter  reaction  used  by  Chapman23  to  describe  the  Na  nightglow.  The  recent 
rate  measurements  of  reaction  (1)  however,  lead  us  to  believe  that  below  80  km 
most  of  the  Na  is  converted  to  NaO  2* 

The  fate  of  Na0  2«  however,  is  uncertain.  Recent  measurements  by  Figger 

24  17 

et  al  and  flame  studies  by  Hynes  et  al  imply  that  the  Na-0  2  bond  strength 
is  much  weaker  than  previously  believed,  with  a  value  of  ca.  163  kJ  mol-1.  As 
a  result,  Na02  will  react  exothermically  with  atomic  oxygen  to  produce  NaO. 
Sodium  oxide  can  react  either  with  ozone  to  reform  Na02  or  react  with  water  to 
produce  NaOH. 

Na02  +  0  +  NaO  +  02  (5) 

NaO  +  03  -►  Na02  +  02  (6) 

NaO  +  H20  -*•  NaOH  +  OH  (7) 

Unfortunately  at  the  beginning  of  this  project,  none  of  these  rates  had  been 
reliably  measured,  leading  the  various  modelers  to  speculate  on  the  importance 
of  each.10*12-15  As  shown  by  Sze  et  al,13  the  predominant  alkali  species 
below  70  km  may  be  NaO,  Na02>  and/or  NaOH,  depending  on  ones  choice  for  kg, 
k6,  and  k7. 

In  a  recent  paper  by  Murad  et  al,8  it  was  proposed  that  the  reactions  of 
metal  hydroxides  and  superoxides  with  chlorine  compounds  between  40  and  70  km 
may  have  an  impact  on  the  depletion  of  stratospheric  ozone.  In  the  case  of 
sodium,  the  exothermic  bimolecular  reactions 

NaOH  +  HC1  +  NaCl  +  ^0  ,  (8) 

NaOH  +  CIO  +  NaCl  +  H02  , 


(9) 


NaOH  +  Cl  NaCl  +  OH  (10) 

Na02  +•  HC1  -►  NaCl  +  H02>  (11) 

Na02  +  CIO  -*•  NaCl  +  03  ,  (12) 

Na02  +  Cl  ->■  NaCl  +  02  (13) 


might  be  expected  to  proceed  rapidly  and  act  as  a  sink  for  Cl,  given  that  NaCl 
can  readily  polymerize  and  condense  via  heterogeneous  nucleation.8  Murad 
et  al  calculated  that  if  the  reaction  rate  of  these  sodium  species  with  HC1 
were  ~10-11  cm3  molecule-1  s-1,  they  might  be  comparable  to  the  major  Cl 
regeneration  mechanism, 

OH  +  HC1  -*•  Cl  +  H20  .  (14) 

While  previously  published  studies  have  viewed  NaCl  as  a  potential  sink  for 

7—8  25 

stratospheric  chlorine,  more  recent  analyses  by  F.S.  Rowland  Indicate 

that  photolysis  of  NaCl  may  in  fact  release  free  Cl.  Given  the  potentially 

large  J  values  (atmospheric  photolysis  rate)  for  this  process, 25-26  reactions 

(8)-(13)  could  effectively  supplement,  rather  than  remove,  reaction  (14)  as  a 

release  mechanism  for  Cl  from  the  Inactive  HC1  stratospheric  reservoir  and 

thereby  determine  the  extent  to  which  ozone  might  be  depleted  by  chlorine 

compounds  in  the  stratosphere.  However,  to  fully  understand  the  role  of 

alkali  species  in  the  stratosphere  one  must  also  consider  the  effects  of  Na02 

and  NaOH  photodissociation  on  those  processes,13’27  as  well  as  homonuclear 

and  heteronuclear  sinks  for  meteoric  metal  species. 

The  objective  of  this  work  is  to  provide  direct  experimental  measurements 
of  these  key  rate  constants.  Obtaining  accurate  reaction  rates  for  the  sodium 
atmospheric  system  will  enable  more  accurate  modeling  of  the  sodium  emission 
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and  density  profiles.  Once  this  is  completed,  these  results  can  be  extended 
toward  the  larger  goal  of  understanding  the  chemistry  of  the  more  abundant 
meteor  metals. 

During  the  first  year  of  this  program,  rate  constants  for  the 
recombination  reactions  of  alkali  atoms  with  molecular  oxygen,  K  +  02+  M  + 
K02  +  M,  and  Na  +  02  +  M  -*■  Na0  2  +  M,  were  measured  as  a  function  of. 
temperature  from  300  to  700  K  using  a  fast  flow  reactor  (see  Appendix  A). 

Laser  induced  fluorescence  was  used  to  monitor  the  disappearance  of  Na  or  K  as 
a  function  of  02  and  M.  The  reactions  were  studied  in  their  low  pressure 
third  order  limit  from  1  to  7  torr  total  pressure  with  N2,  He,  and  Ar  as  third 
bodies,  and  the  rate  constants  were  found  to  have  the  expected  negative 
temperature  dependence.  These  results  compared  well  with  other  recent 
measurements  from  flame  and  flash  photolysis  studies  and  with  theoretical 
expectations  based  on  an  energy  transfer  RSKM  mechanism  for  the  Na02* 
activated  complex.  Similar  techniques  were  used  to  demonstrate  that  the 
reactions  of  NaO  and  NaOH  with  HC1  had  a  gas  kinetic  rate  constant.  (Appendix 
B)  In  addition,  limits  for  the  rate  constants  were  found  for  the  reactions  of 
atomic  hydrogen  with  NaO,  NaOH,  and  NaCl. 

The  work  during  this  past  year  has  continued  and  expanded  our  efforts  of 
the  first  year.  We  have  measured  reaction  rate  constants  for  Na  +  O3  and  NaO 
+  03  at  294  K,  and  have  preliminary  data  for  Na02  +  HC1  and  for 
photodissociation  cross  sections  of  NaCl.  In  addition,  we  have  begun  to 
examine  the  various  ways  in  which  the  reaction  neutral  gas  phase  alkali 
species  are  removed  from  the  mesosphere  and  stratosphere. 

1 .2  Work  Forecast 

During  the  next  year  of  this  program,  we  will  continue  our  examination  of 
the  chemistry  of  alkali  metals  in  the  mesosphere  and  upper  stratosphere.  The 
major  reactions  of  interest  are  NaO 2  +  HC1,  NaO  +  H2O  and  NaOH  +  CO  2  +  M.  In 
addition,  we  hope  to  measure  photodissociation  cross  sections  for  NaOH,  NaCl, 
NaO ,  and  NaO  2  • 
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2 .  TECHNICAL  PROGRESS 


I 


2.1  Introduction 

We  have  made  good  progress  during  this  second  year  of  our  program  to 
study  the  gas  phase  kinetics  of  meteoric  metals  in  the  atmosphere.  In  this 
section,  we  will  describe  in  detail  the  experimental  results  obtained  this 
year  and  mention  our  plans  for  the  coming  year.  The  experiments  have  been 
performed  in  the  Aerodyne  high  temperature,  fast  flow  reactor,  which  has  been 
described  in  detail  previously 18» 28-29  and  will  not  be  repeated  here. 


2 .2  Measurements  of  Atomic  Sodium  and  Potassium  Diffusion  Coefficients 

In  the  process  of  analyzing  experimental  kinetic  data,  it  is  necessary  to 
make  corrections  for  diffusion  in  the  observed  first-order  decay  rates.  In 
our  experiments,  the  probability  of  any  alkali  species  sticking  to  the  wall  of 
the  flow  tube  is  unity.  This  fact  allowed  us  (see  Appendix  C  for  details)  to 
directly  determine  diffusion  coefficients  from  observed  wall  loss 
measurements,  and  resulted  in  a  publication  in  the  Journal  of  Chemical 
Physics. 90 

The  results  are  shown  in  Table  2.1. 


Table  2.1  -  Diffusion  Coefficients  for  Na  and  K  in  Various  Gases 


GASES  D( cm2  s-1  at  1  atm) 

Na-He  (2. 2+0. 3)  x  10-5  x1’75!0’02 

Na-N2  (8. 0+0. 8)  x  10“ 6  x1’79*0*02 

Na-Ar  (1.2+0. 9)  x  10-5  T1’71*0'  11 

K-He  (1. 5+1.0)  x  10" 3  T1*^0*  19 

K-N2  (4. 5+2. 8)  x  10-5  x1*54*!0’09 


T  (K) 
309-473 
320-698 
322-350 

301- 734 

302- 720 


7 


2.3  Reaction  Rate  Constants  for  Na  and  NaO  with  0 3 


2.3.1  Ozone  Production  and  Detection 

Measurement  of  the  reaction  rate  constant  of  atomic  sodium  with  ozone  is 
straightforward,  given  a  quantifiable  source  of  ozone.  The  loss  rate  of 
atomic  sodium  is  monitored  by  laser  induced  fluorescence  as  discussed 
previously. 29  In  these  studies,  03  is  in  excess  of  Na  so  that  a  pseudo-first 
order  kinetic  analysis  can  be  used. 

Ozone  is  produced  in  a  small  commercial  ozonator  (Uelsbach  Model  T-408) 
and  stored  by  cryogenic  trapping  on  silica  gel  in  a  specially  designed  trap. 
Once  stored  in  this  trap,  ozone  is  Introduced  into  the  flow  tube  by  eluting  it 
with  a  known  flow  of  helium.  The  concentration  of  0 3  in  this  flow  is 
determined  using  a  compact,  calibrated  absorption  cell,  similar  in  design  to 
that  described  in  Ref.  31.  The  partial  pressure  of  03  is  calculated  using  the 
measured  cross  section  for  the  coincidental  absorption  of  the  253.7  nm  mercury 
line,  generated  in  a  standard  pen  lamp.  The  flow  rate  of  ozone  is  then 


3  -1 

(STP  cm  s 


) 


( 


Pto«l  -  V  <STP  C"V1) 


(15) 


2.3.2  Rate  Measurements  and  Analysis 

In  preparing  to  measure  the  rate  constant  of  the  reaction 


Na  +  03  -*■  NaO  +  O2  AH°  *  -39.9  kcal/mole,  (2) 

it  turns  out  that  we  must  also  be  aware  of  possible  complications  due  to 
secondary  reactions  of  the  NaO  product, 


28.4  kcal/mole,  and  (4a) 


NaO  +  03  +  Na  +  20z 


NaO  +  O3  -*•  NaO 2  +  02 


AH°  -  -  67.4  kcal/mole.  (4b) 


The  termolecular  combination  reaction, 

Na  +  02  +  M  -*•  NaO 2  +  M,  (1) 

is  too  slow2  at  the  experimental  pressures  to  be  of  importance.  If  the 
products  of  reaction  (4)  are  Na02  +  02,  then  no  complications  In  analyzing  the 
[Na]  decay  for  k2  would  exist.  However,  if  reaction  (4a)  plays  a  significant 
role,  the  Na  formed  in  this  reaction  will  subsequently  react  with  03,  causing 
the  observed  [Na]-time  profiles  to  depart  from  a  simple  logarithmic  decay.  An 
example  of  predicted  behavior  for  [Na]  decays  is  illustrated  in  Figure  2-1. 

In  this  calculation,  assume  that  k2  -  4  x  10“ 10  cm3s“1  and  vary  both  the 
value  of  k4  and  the  ratio  k4a/k4ij.  As  expected,  if  k4a  ■  0,  then  the 
observed  decay  corresponds  to  k2.  If  k4a  is  finite,  then  the  various  curves 
in  [Na]  will  be  observed. 

The  mathematical  expression  which  describes  the  curves  in  Figure  1  can  be 
written  as 


(k4-k2+A)e 


-[03] (k2+k4-A) t/2 


-  (k4-k2-A)e 


-[°3](k2+k4+A)t/2 


[Na( t) ] 


2A 


(16) 


where  k4  -  k4a  +  k4b  and 
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4  x  10 


5  x  10 


2.5  x  10 


5  x  10 


2.5  x  10 


2.5  x  10 


2.5  x  10 


2.5  x  10 


2.5  x  10 


Figure  2-1.  Calculated  Sodium  Profiles  for^Var^ous^Values  of  k^a 
and  k^*  Assuming  =  4  x  10  cm  s 


^  m  [(^2  +  )2  ”  4k-2  ^14. b] 


1/2 


Thus,  the  Na  concentration  varies  as  the  sum  of  two  exponential  terms  whose 
variables  are  [03]t.  We  can  then  combine  all  of  our  data,  plotting  the  points 
as  log[ Na]  vs.  [03]t,  and  obtain  estimates  of  k4a  and  k4^* 

2.3.3  Results 

A  series  of  measurements  were  performed,  and  curved  decay  plots  were 
observed.  Indicating  that  the  value  of  k1+a  is  not  negligible.  In  order  to 
more  accurately  obtain  kj  ,  we  added  an  excess  of  HC1  to  the  flow.  This  has 
the  effect  of  suppreslng  reactions  4a  and  4b.  The  reaction 

NaO  +  HC1  ♦  NaCl  +  OH  (17) 

will  overwhelm  reaction  4  because  -  2.8  x  10-10  cm3  molecule-1  s-1  and 
[HC1J  »  [03 ] .  We  observe  very  linear  decays  over  many  orders  of  magnitude 
in  signal  (Fig.  2-2).  Correction  of  these  data  for  diffuslonal  and  wall 
effects  leads  to  a  value  of  1*2  ■  3.2  +  0.4  x  10-10  cm3  molecule-1  s-1 ,  where 
the  uncertainty  is  one  standard  deviation  and  includes  consideration  of  both 
random  and  systematic  errors  (Fig.  2-3).  A  simple  prediction  for  this 
reaction  rate,  based  on  an  electron  jump  model,33  would  predict  a  value  of  4.7 
x  10-10  cm3  molecule-1  s-1  .  A  recent  measurement31*  using  a  similar  approach 
as  ours  obtained  a  value  for  k2  of  7  x  10-10  cm3  molecule-1  s-1.  We  are  in 
the  process  of  trying  to  resolve  the  factor  of  two  difference  between  these 
results. 

Substituting  k2  into  equation  (16),  we  then  attempted  to  fit  our  data  to 
obtain  values  for  lq,a  and  k^*  Unfortunately,  the  scatter  in  the  data  is 
too  large  at  this  time  to  get  answers  to  better  than  a  factor  of  two.  As  the 
value  of  k^a  +  k^  approaches  or  exceeds  that  of  k4 ,  the  fit  becomes 
insensitive  to  value  of  kjg  +  k^i  especially  if  k^^  >  k4a.  However, 
based  on  the  overall  shape  of  the  curves,  we  can  deduce  that  k4a  ■  3k4^ 
and  that  the  overall  reaction  rate  constant  for  NaO  +  03  is  >10-10  cm3 
molecule-1  s-1 . 
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2.4  Rate  Constant  for  NaO?  +  HC1 


Ve  have  begun  measurements  to  determine  the  rate  constant  for  the 
reaction 


Na02  +  HC1  -*•  NaCl  +  H02 


Sodium  superoxide  is  formed  in  the  upstream  region  of  the  flow  tube  by  the 
termolecular  reaction 

Na  +  02  “  Na02,  *  (1) 

and  is  detected  in  a  manner  described  in  Appendix  B,  whereby  the  Na02  is 
titrated  in  the  detector  region  by  atomic  hydrogen,  and  the  resulting  sodlumn 
atom  detected  by  laser-induced  fluorescence. 


Na02  +  H  -►  Na  +  H02 


Na  +  hv  ♦  Na  Na  +  hv 


(LIF) 


Because  the  NaCl  product  also  forms  Na  upon  addition  of  H, 
NaCl  +  H  -*•  Na  +  HC1, 


(albeit  at  a  slower  rate),  the  observed  decay  plots  are  curved  and  eventually 
level  off  at  a  value  reflecting  the  relative  rate  constants  of  reactions  18 
and  20. 

Preliminary  results  Indicate  that  the  reaction  of  Na02  with  HC1  is 
relatively  fast.  Although  more  experiments  are  needed,  at  this  time  a  lower 
limit  of  k  can  be  set  at  3  x  10” 11  cm3  molecules"1  s"1  at  294  K. 

2.5  Photodlssoclatlon  Cross  Sections  for  NaCl 

2.5.1  Introduction 

If  alkali  chemistry  is  going  to  have  an  impact  on  ozone  reduction  in  the 
upper  stratosphere,  three  criteria  must  be  met:  1)  Na  must  be  rapidly 


transformed  to  NaCI,  2)  photolysis  must  quickly  regenerate  free  Cl  and  atomic 
Na  from  NaCI,  and  3)  homogeneous  and  heterogeneous  processes  which  remove  the 
mononeric  neutral  gas  phase  alkali  molecules  must  not  be  too  rapid.  It 
appears  that  the  first  criteria  is  met,  and  we  will  discuss  the  third  in  the 
next  section  of  this  report.  The  photodissociation  cross  sections  of  alkali 
molecules,  and  most  importantly  NaCI,  need  to  be  measured  accurately  between 
approximately  200  nm  and  300  nm,  the  latter  value  corresponds  to  just  below 
the  dissociation  band  energy  of  NaCI. 

Although  a  few  previous  measurements  of  alkali  halide  photoabsorption 
cross  sections  have  been  made,35-37  they  all  (with  one  exception37  for  Csl), 
have  been  done  at  higher  temperatures,  where  dimers  and/or  vibrational  hot 
bands  may  have  significantly  contributed  to  the  observed  spectra.  Since  the 
photolysis  rate  is  a  convolution  of  the  solar  spectrum  (as  seen  in  the  upper 
stratosphere)  and  the  photodissociation  cross  section,  it  is  desireable  to 
have  detailed  cross  sections  which  are  free  of  the  effect  of  dimers  and 
correspond  to  a  stratospheric  vibrational  temperature  of  250  K.  Our  approach 
of  making  NaCI  at  room  temperature  by  the  chemical  reaction  of  Na  +  Cl2  more 
closely  satisfies  these  requirements  and  should  result  in  an  accurate 
determination  of  the  stratospheric  photolysis  rate  of  NaCI. 


2.5.2  Experimental  Approach  and  Preliminary  Results 

Sodium  chloride  is  formed  by  the  rapid  reaction  of  Cl  2  +  Na.  We  have 
measured  this  rate  constant  to  be  '6  x  10" 10  cm3  molecule" 3  s" *  at  294  K. 

Our  initial  approach  toward  measuring  the  cross  section  was  to  use  a  high 
Intensity  Xe  lamp,  dispersed  through  a  monochromator,  as  our  light  source.  As 
NaCI  flowed  through  an  irradiated  region,  it  would  photodissociate  and  the 
amount  of  detected  Na  just  downstream  relative  to  the  initial  amount  of  NaCI 
will  be  expressed  as 


InScT]o 


l-e-oUUoUH 


(21) 
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where  IQ  -  photons  cm  s“  at  wavelength  \ 
a(\)  m  cross  section  in  cm2  at  \ 

t  *  Irradiation  time  in  seconds  *  iradiation  length  (d)/gas 

velocity  (v) 

Thus  a  measurement  of  the  left-hand  side  of  equation  21,  coupled  with  an 
accurate  measurement  of  10»  d*  an^  v»  would  result  in  <7  (\). 

Several  problems  arose,  however,  in  performing  these  measurements. 

First,  the  geometry  of  the  flow  tube  is  such  that  the  flow  velocity  in  the 
detection  region  is  not  equal  to  that  in  the  reaction  portion.  Although  this 
has  no  effect  on  kinetic  measurements,  it  means  that  the  flow  velocity  in  this 
region  is  not  easily  calculated,  but  must  be  measured.  Although  possible  by  a 
variety  of  techniques,  it  requires  additional  work.  The  second,  more  serious 
problem  was  the  Xe  illuminator.  The  500  Watt  lamp  we  had  was  incapable  of 
being  focussed  properly  through  the  monochromator  without  losing  too  much 
intensity.  Although  a  smaller  75  Watt  lamp  of  the  proper  design  was  used,  it 
was  too  weak  to  provide  measureable  dissociation  cross  sections. 

After  considering  the  alternate  choice  of  acquiring  a  new  1000  Watt  lamp, 
we  decided  to  instead  go  to  a  laser  photolysis  setup.  In  this  arrangement,  a 
Raman  shifted  excimer  laser  will  replace  the  Xenon  lamp.  This  approach  has 
several  advantages.  The  laser  pulse  is  only  about  10  ns  long,  so  that  the 
NaCl  molecules  are  effectively  stationary  during  photolysis.  Thus  only  the 
total  photon  density  in  each  pulse  is  required  to  be  known  (i.e.,  pulse  energy 
and  laser  beam  diameter);  we  no  longer  need  to  know  v.  Secondly, 
instantaneous  pulse  energy  is  higher  at  any  selected  wavelength,  so  that 
[Na Jj/ [NaCl ]Q  will  be  larger  and  more  easily  measured  ( ~10“ 2  -  10“ 3 
rather  than  lO-*4  -  10"6  with  the  500  W  Xe  lamp).  The  probe  visible  laser  is 
now  coincident  with  the  pump  (UV)  laser  and  detects  the  photodissociated  Na 
within  a  microsecond  of  formation.  NaCl  is  measured  by  turning  the  Cl2  and 
excimer  laser  off,  making  the  assumption  that  all  of  the  original  sodium  is 


completly  converted  to  NaCI  upon  addition  of  excess  Cl  2*  Since  the  experiment 
is  over  in  <1  ps,  there  is  insufficient  time  for  the  newly  formed  sodium  atoms 
to  react  with  this  remaining  Cl 2  before  they  are  detected.  Using  H2/D2 
mixtures  in  the  Raman  cell,  we  can  easily  cover  the  200-300  nm  region  with  at 
least  25  wavelengths . 

At  this  time  we  have  demonstrated  the  viability  of  this  approach  using 

the  excimer  laser  alone  at  193  nm.  Preliminary  analysis  gives  0=  (1.3+0. 6)  x 
17  2 

10"  cm  at  this  wavelength.  Preparations  are  now  underway  for  extending 
these  measurements  to  300  nm  using  the  Raman  shifter. 

2.6  Atmospheric  Sinks  for  Meteoric  Alkali  Compounds 

As  demonstrated  above,  alkali  gas  phase  chemistry  may  significantly 
impact  chlorine  induced  ozone  depletion  if  0.1-1%  of  the  potential 
stratospheric  alkali  content  remains  in  monomeric  molecular  form. 

Various  sinks  for  gaseous  alkalis  can  be  considered.  Let  us  look  at 
three  potential  sinks  for  neutral  alkali  species: 

2.6.1  Ionization/Cluster  Formation 

Hydration  of  alkali  species  by  H30T*"»nH2O  ions  will  remove  neutral  alkali 
and  form  ions  which  rapidly  cluster  and/or  condense  on  aerosols. 


NaX  +  H30*-»nH20  ->•  NaXH*-  »nH  20  +  H20 

(22) 

NaXH*-.nH20  +  NaX  -*•  (NaX)  2.(n-l)H20  +  H20 

(23) 

[NaX][H20]  >  1/K  (23) 

(24) 

Upper  stratospheric  balloon-borne  mass  spectrometer  measurements  by  Arnold  et 
al38and  Arijs  et  al39  show  no  identifiable  NaXH^nH^O  peaks,  thus  setting 
limits  on  the  amount  of  neutral  alkali  present.  Subject  to  better  refinements 
of  Kgq  and  mass  spectrometer  sensitivities,  this  present  limit  is  not 
incompatible  with  an  NaX  of  -10 **  cm-3. 148  Furthermore,  the  photolytic 
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stability  of  either  Ion  or  neutral  clusters  to  reform  NaX  is  unknown.  Recent 
unpublished  results  from  rocket  launched,  parachute  deployed  mass  spectroscopy 
measurements  of  upper  stratospheric  and  lower  mesospheric  ion  measurements  by 
Arnold  do  indicate  a  strong  possibility  of  either  sodium  or  magnesium  core 
ions  in  the  40  to  60  km  range.  ^ 

2.6.2  Homogeneous  Condensation 

Dimerization  and  subsequent  clustering  of  alkali  species  will  remove  them 
from  the  catalytic  Cl  regeneration  cycle  discussed  above.  At  night,  the 
alkalis  will  probably  reside  in  the  form  of  NaCl,  which  in  the  absence  of 
light  will  dimerize. 

NaCl  +  NaCl  (NaCl)2 

Even  with  an  extremely  fast  termolecular  recombination  rate  of  10“ 28  cm6 
2  1  3  3 

molecule”  s-  ,  more  than  5  x  10  cm”  NaCl  will  still  be  in  the  monomeric 
form  after  4  months,  which  is  the  upper  limit  of  estimated  replacement  time  of 
sodium  in  the  stratosphere  by  diffusion  (10 6  -  10 7  s).1*2 

2.6.3  Heterogeneous  Condensation 

It  has  been  suggested  by  Hunten  and  co-workers 42  that  ablation  of  meteors 
may  lead  the  formation  of  high  altitude  "smoke",  small  particles  with  an 
integrated  equivalent  surface  area  (a)  of  10”9  cm2/cm3  between  30  and  85  km. 

As  the  alkali  species  diffuse  downward,  they  might  condense  on  these  aerosols 
and  be  removed  from  the  gas  phase.  The  amount  of  loss  depends  critically  on 
the  accommodation  (sticking)  coefficient,  a,  and  on  the  residence  time  in  the 
atmosphere  between  85  and  30  km. 

—  /VflWf 

(NaX)  -  (NaX)Qe 

With  gas-particle  relative  velocities  (v)  of  ~3  x  10  **  cm  s”1,  [NaX  ]/ [NaX  ]0 
ranges  from  0-0.7  for  a  in  the  range  1.0-0.01  and  t  *  106-107  s.  Not  only  are 
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the  accommodation  coefficients  unknown,  but  the  net  loss  also  depends  on 
whether  surface  photolysis  and/or  chemistry  releases  adsorbed  NaX  species. 

2.7  Summary 

In  the  first  two  years  of  this  program,  much  has  been  learned  about  the 
processes  governing  the  chemistry  of  meteoric  metal  in  the  atmosphere.  The 
chemical  fate  of  Na  in  the  lower  thermosphere  is  being  clarified  due  to  our 
results  on  Na  +  O2  +  M,  Na  +  O3,  and  NaO  +  O3.  Further  work  on  NaO  +  H2O  and 
NaOH  +  CO 2  +  M,  as  well  as  loss  mechanisms,  should  provide  a  firm  basis  for 
predicting  species  profiles.  The  stratospheric  chemistry  of  meteor  metals  has 
been  better  understood  due  to  our  results  for  NaO,  NaOH,  and  Na02  with  HC1. 
Further  work  in  this  area  will  involve  the  determination  of  photodissociation 
cross  sections  for  the  alkali  oxides,  hydroxides,  and  halides. 
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6.  PRESENTATIONS 


MEETINGS  AND  TRAVEL 


During  this  past  year,  the  following  presentations  were  made: 


"The  Reaction  Rates  of  Na  +  Oj  +  M,  K  +  O2  +  M  and  NaOH  +  HC1”,  by  J.A. 
Silver,  M.S.  Zahniser,  A.C.  Stanton,  and  C.E.  Kolb,  at  the  1983  Fall 
American  Geophysical  Union  meeting  in  San  Francisco,  CA  on  December  7, 

1983. 

"Atmospheric  Chemistry  of  Meteoric  Metals",  by  J.A.  Silver  for  a 
colloquium  at  the  Air  Force  Geophysics  Laboratory  on  April  19,  1984. 

"Temperature  Dependent  Termolecular  Reaction  Rate  Constants  for  Potassium 
and  Sodium  Superoxide  Formation",  by  J.A.  Silver,  M.S.  Zahniser,  A.C. 
Stanton,  and  C.E.  Kolb,  at  the  20th  International  Symposium  on  Combustion 
in  Ann  Arbor,  MI  in  July,  1984. 

"The  Chemistry  of  Alkali  Oxides  in  the  Atmosphere",  by  J.A.  Silver  at  the 
XVI  Informal  Conference  on  Photochemistry  in  Cambridge,  MA  in  August, 

1984. 

"Gas  Phase  Alkali  Chemistry  in  the  Post-Molecular  Beam  Age",  by  C.E.  Kolb 
at  MIT  on  6  November  1984. 
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and  Sodium  Superoxide  Formation",  by  J.A.  Sliver,  M.S.  Zahniser,  A.C. 
Stanton,  and  C.E.  Kolb,  20th  Symposium  (Int.)  on  Combustion,  The 
Combustion  Institute  (1984). 

"The  Gas  Phase  Reaction  Rate  of  Sodium  Hydroxide  with  Hydrochloric  Acid”, 
J.A.  Silver,  A.C.  Stanton,  M.S.  Zahniser,  and  C.E.  Kolb,  J.  Phys.  Chem., 
88,  3123  (1984). 

"Measurement  of  Atomic  Sodium  and  Potassium  Diffusion  Coefficients", 

J.A.  Silver,  J.  Chem.  Phys.,  81,  5125  (1984). 


26 


APPENDIX  A 


TEMPERATURE  DEPENDENT  TERMOLECULAR  REACTION  RATE  CONSTANTS  FOR  POTASSIUM 

AND  SODIUM  SUPEROXIDE  FORMATION 


To  be  published,  20th  Symposium  (International)  on  Combustion,  1984 


1 


ARI-RR-381 


TEMPERATURE  DEPENDENT  TERMOLECULAR  REACTION 
RATE  CONSTANTS  FOR  POTASSIUM  AND 
SODIUM  SUPEROXIDE  FORMATION 


Prepared  by 


J.A.  Silver,  M.S.  Zahniser,  A.C.  Stanton,  and  C.E.  Kolb 
Center  for  Chemical  and  Environmental  Physics 
Aerodyne  Research,  Inc. 

45  Manning  Road 
Billerica,  MA  01821 


Prepared  for 


The  20th  Symposium  (International)  on  Combustion 
The  Combustion  Institute 
Pittsburgh,  PA 


January  1984 


Subject  Matter:  Kinetics 


28 


ABSTRACT 


Race  constants  for  the  recombination  reactions  of  alkali  atoms  with 
molecular  oxygen,  K  +  02  +  M  +  K02  +  M  and  Na  +  02  +  M  ■*  Na02  +  M,  have  been 
measured  as  a  function  of  temperature  from  300  to  700  K  using  a  fast  flow 
reactor.  Laser  induced  fluorescence  is  used  to  monitor  the  disappearance  of 
Na  or  K  as  a  function  of  02  and  M.  The  reactions  are  studied  in  their  low 
pressure  third  order  limit  from  1  to  8  Corr  total  pressure  with  N2,  He,  and  Ar 
as  third  bodies.  The  rate  constants  are  found  to  have  the  expected  negative 
temperature  dependence.  The  values  for  k  (Na  +  02  +  M)  are  (1.9  ±0.4)  x 
10"30  (T/300)”1 * 1  ±0*5  with  M  -  N2,  (1.4  ±  0.3)  x  10-30  (T/300)-0*9  ±0*5 
with  M  -  He,  and  (1.2  ±  0.3)  x  10-30  at  T  ■  324  K  with  M  -  Ar,  all  in  units  of 
cm6  molecules”2  s-1.  The  values  for  the  corresponding  reactions  with 
potassium  are  larger  than  those  for  sodium  with  k  (K  +  02  +  M)  -  (5.4  ±  0.2)  x 
10-3°  (x/300)“°-56  ±0*20  with  M  -  N2,  (2.0  ±  0.5)  x  10-30  (T/300)-0*9  ±°*5 
with  M  ■  He,  and  (3.5  ±1.0)  x  10-30  at  T  ■  300  K  with  M  -  Ar.  The  results  are 
compared  with  other  recent  measurements  from  flame  and  flash  photolysis 
studies  and  with  theoretical  expectations  based  on  an  energy  transfer  RRKM 
mechanism  for  the  Na02*  activated  complex. 
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INTRODUCTION 


The  chemistry  of  alkali  metals  in  flames  commands  our  attention  for  both 
scholarly  and  practical  reasons.  From  a  scientific  standpoint  reactions  of 
alkali  atoms,  with  their  single  valence  electron,  form  a  natural  test  bed  to 
extend  our  theoretical  understanding  honed  on  reactions  of  atomic  hydrogen. 
Furthermore,  the  low  ionization  potentials  of  alkali  atoms  open  the 
possibility  of  electron  transfer  reaction  mechanisms  and  allow  determination 
of  the  role  of  ionic  potential  surfaces  in  reaction  dynamics.  Finally, 
incredibly  sensitive  methods  of  detecting  alkali  atoms,  including  hot  wire 
surface  ionization,  laser  induced  resonance  fluorescence,  and  laser  induced 
resonance  ionization  now  allow  the  design  of  elegant  experimental  studies  for 
environments  ranging  from  molecular  beams  to  high  pressure  flames. 

From  a  practical  standpoint,  gas  phase  alkali  chemistry  has  long  been 
recognized  to  be  important  in  flame  suppression, 1-5  and  is  also  involved  in 
the  conversion  of  alkali-containing  minerals  in  coal  to  alkali  sulfates.6-8 
These  sulfates  are  a  major  cause  of  fouling  and  corrosion  of  boiler-tube 
surfaces,  heat  exchangers,  and  turbine  blades.  Alkali  atoms  introduced  into 
the  atmosphere  by  meteor  ablation  also  play  a  role  in  mesospheric 
chemistry,9-12  and  it  has  been  suggested13-11*  that  alkali  molecules  may  affect 
stratospheric  ozone  reduction  through  the  catalytic  chlorine  cycle. 

However,  for  all  the  intrinsic  theoretical  and  experimental  interest  in 
alkali  reactions,  our  understanding  of  the  gas  phase  oxidation  chemistry  of 
alkali  atoms  has  only  recently  progressed,  first  back  to,  and  then  beyond  the 
level  gained  in  the  1930's  in  pioneering  sodium  diffusion  flame  studies  by 
Haber  and  Sachsse15  and  Bawn  and  Evans.16  While  these  low  pressure  studies  of 
reaction  (2)  showed  it  proceeded  with  a  termolecular  rate  constant  in  excess 
of  10-3°  cm6  s-1,  subsequent  work  (primarily  in  the  mid  1960’s)  using  flame 
photometric  methods  by  Kaskan17-19  and  McEwan  and  Phillips 19-20  indicated  much 
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lower  M02  formation  rates  (M-Na,K) 17-19  as  well  as  relatively  high 
dissociation  energies  for  the  M-02  bond  (M«Na, Li). 19-20 

Analyses  of  thermochemical  cycles  by  Herm  and  Herschbach21  and 
Alexander22  clearly  indicated  that  the  values  of  tt-02  bonds  deduced  from  flame 
photometric  data  by  McEwan  and  Phillips  were  greatly  overestimated.  This  was 
recently  confirmed  by  Flgger  and  co-workers  who  analyzed  observations  of 
chemiluminescence  between  crossed  molecular  beams  of  alkali  dimers  and 
molecular  oxygen  to  yield  alkali  superoxide  bond  energies.23  Reinterpretation 
of  the  earlier  flame  experiments2**  as  well  as  newly  performed  flame  studies23 
are  also  consistent  with  lower  alkali  superoxide  bond  energies. 

The  low  sodium  and  potassium  termolecular  reaction  rates  deduced  from  the 
flame  photometry  studies17-19  are  also  incorrect.  Recently  published  flash 
photolysis  studies  over  limited  temperature  ranges  for  Li  (393  and  463  K),26 
Cs  (323  K),27  Na  (724  and  844  K),28  and  K  (753  and  853  K),29  as  well  as  flame 
photometry  studies  for  Na,25  all  clearly  show  fast  termolecular  rates  for 
alkali  superoxide  formation. 

In  this  paper  we  report  termolecular  rate  constants  for  superoxide 
formation  of  sodium  and  potassium  over  a  much  wider  temperature  range  than 
previously  available,  yielding  the  first  clear,  direct  measurement  of  the 
temperature  dependence  of  these  processes.  Third  body  dependencies  are 
illustrated  by  separate  measurements  for  N2,  He  and  Ar.  These  measurements 
were  made  in  a  high  temperature,  fast  flow  reactor30  utilizing  laser  induced 
resonance  fluorescence  detection  of  atomic  alkali  species.  A  discussion  of 
these  results  in  terms  of  their  importance  in  testing  termolecular  association 
reaction  theory  is  also  presented. 
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EXPERIMENTAL 


The  high  temperature  fast  flow  reactor  used  in  these  experiments  has  been 
described  in  detail  previously30  and  is  shown  in  Fig.  1  in  the  configuration 
used  in  this  study.  Briefly,  a  7.26  cm  diameter,  120  cm  long  alumina  tube  is 
heated  with  Kanthal  resistance  heaters.  The  helium,  argon  or  nitrogen  carrier 
gas  enters  the  flow  tube  through  two  mulllte  multichannel  arrays  which  both 
preheat  and  laminarize  the  flow  upstream  of  the  reaction  zone.  A  sufficient 
distance  (20  cm)  is  allowed  downstream  of  these  arrays  for  the  carrier  gas  to 
develop  a  parabolic  flow  profile  before  reaching  the  reaction  zone.  Carrier 
gas  flow  rates  are  measured  with  calibrated  rotameters.  Total  pressures  are 
measured  with  a  capacitance  manometer  at  the  downstream  end  of  the  reaction  . 
zone.  Gas  temperatures  are  obtained  with  a  chromel-alumel  thermocouple  which 
is  movable  throughout  the  reaction  zone.  The  lowest  reaction  temperatures  are 
slightly  above  room  temperature  due  to  the  heating  of  the  carrier  gas  as  it 
passes  over  the  resistively  heated  alkali  oven.  The  maximum  axial  temperature 
gradient  due  to  this  effect  over  the  region  in  which  decay  measurements  are 
taken  is  10  to  20  degrees  depending  on  flow  velocity,  pressure,  and  identity 
of  carrier  gas.  The  temperature  at  the  mid  point  of  the  reaction  zone  as 
measured  by  the  movable  thermocouple  is  used  In  the  data  analysis.  Axial 
profiles  are  more  uniform,  ±2  to  5  degrees,  at  elevated  temperatures  (>  400  K) 
where  heating  is  dominated  by  the  flow  tube  walls  and  the  mulllte  arrays. 
Although  the  apparatus  is  capable  of  temperatures  up  to  1500  K,  experimental 
considerations  as  described  below  limited  the  upper  temperature  to  700  K  for 
these  studies. 


Alkali  atoms  are  generated  by  heating  the  pure  metal  in  a  25  mm  diameter 
cylindrical  silver  plated  monel  oven  to  a  temperature  sufficient  to  obtain  a 
vapor  pressure  of  10"6  to  10”1*  torr  within  the  oven.  The  oven  is  mounted  on 
the  end  of  a  movable  13  mm  od  alumina  tube  concentric  to  the  main  flow  tube. 
The  vapor  is  entrained  in  a  flow  of  carrier  gas  and  introduced  into  the  flow 
tube  in  one  of  two  configurations.  For  the  kinetic  measurements  with  02>  the 
oven  is  placed  at  a  fixed  position  upstream  of  the  multichannel  arrays  and  the 
alkali  vapor  is  introduced  into  the  main  flow  through  a  10  cm  length  of  19  mm 
diameter  silver  tubing  which  passes  through  the  center  of  the  arrays.  In  the 
second  configuration,  the  oven  is  placed  downstream  of  the  arrays  and  is 
movable  throughout  the  reaction  zone  to  determine  the  loss  of  alkali  atoms  on 
the  reactor  walls.  The  oven  temperature  is  controlled  with  resistive  heating 
elements  Independent  of  the  flow  tube  heaters.  Initial  alkali  atom 
concentrations  in  the  reaction  zone  are  maintained  at  less  than  10 10  cm~3  by 
adjusting  the  oven  temperature  and  the  carrier  gas  flow  rate  through  the  oven. 

The  alkali  atoms  are  detected  by  laser  Induced  fluorescence  at  the 
downstream  end  of  the  flow  tube.  A  Molectron  DL14  nitrogen  pumped  dye  laser 
is  used  to  excite  the  4s  2$i/2  *  5p  2p3/2  transition  at  404.4  nm  for 
potassium  or  the  3s  2Si/2  2^3/2  transition  at  589.0  nm  for  sodium. 

Fluorescence  is  collected  using  a  gated  phototube  and  a  computer  controlled 
data  acquisition  system.30  The  signals  are  averaged  over  100  laser  pulses 
after  subtracting  for  nonfluorescent  background  contributions  and  normalizing 
for  pulse  to  pulse  fluctuations  in  laser  intensity.  Although  no  direct 
calibration  of  the  fluorescence  signal  was  attempted,  estimates  of  sensitivity 
from  known  phototube  and  integrator  response,  measured  laser  power,  and  atomic 
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transition  probabilities  indicate  a  detection  limit  for  Na  of  104  cm“3  and  for 
K  of  106  cm"3  using  these  transitions. 

Molecular  oxygen  is  added  through  a  movable  loop  injector  at  distances 
from  10  to  60  cm  upstream  of  the  detection  region.  02  flow  rates  are  measured 
with  a  calibrated  thermal  conductivity  mass  flow  meter.  The  kinetic 
measurements  are  conducted  under  pseudo-first  order  conditions  with  2  x  1014 
<[02]<  7  x  1015  cm-3,  in  large  excess  compared  to  the  alkali  metal  atom 
concentration,  yet  always  <5Z  of  the  total  gas  concentration.  First  order 
rate  constants  are  determined  for  a  fixed  02  concentration  by  recording  the 
change  in  alkali  atom  concentration  as  a  function  of  02  injector  distance. 

Corrections  for  both  axial  and  radial  diffusion  and  wall  removal  are  made 
using  the  procedure  outlined  by  Brown.31  This  method  is  based  on  a  numerical 
solution  of  the  continuity  equation  including  diffusion,  first  order  chemical 
reaction  and  wall  removal  in  a  cylindrical  tube  with  fully  developed  laminar 
flow.  A  multiplicative  correction  factor,  obtained  for  each  set  of 
experimental  conditions,  is  used  to  relate  the  observed  decay  to  the  true 
first  order  rate  constant.  This  correction  factor  is  dependent  on  the  wall 
removal  rate  constant,  1^,  and  Che  binary  diffusion  coefficient  for  the 
alkali  atom  in  the  carrier  gas.  Values  for  ky  are  determined  in  separate 
experiments  at  each  pressure  and  temperature  by  observing  the  change  in  atom 
concentration  at  the  detection  point  while  varying  the  source  oven  position. 
The  loss  of  atoms  to  the  wall  was  found  to  be  diffusion  limited  for  all  the 
conditions  used  in  this  study,  implying  a  surface  accommodation  coefficient  y 
>  0.1.  The  observed  wall  loss  rate  under  these  conditions  may  therefore  be 
used  to  determine  the  diffusion  coefficient 3 ^  which  are  then  used  to  determine 
the  correction  factors  for  the  reactive  flow  analysis.  The  correction  factors 
obtained  by  this  method  were  in  the  range  ktrue/k0^g  from  1.52  to  1.86, 
with  moat  lying  between  1.60  and  1.65.  Under  the  flow  conditions  of  these 
experiments,  the  ratio  of  reaction  to  diffusion  times  is  not  unlike  those 
encountered  in  ion-molecule  reactions,33  which  also  exhibit  correction  factors 
on  the  order  of  1.6. 
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RESULTS  AND  DISCUSSION 


Reaction  rate  constants  of  atomic  sodium  with  0 2  were  measured  over  the 
temperature  range  of  320  -  698  K  for  nitrogen  as  the  third  body  collision 
partner,  and  over  the  range  of  309  -  473  K  In  helium;  for  potassium,  the 
corresponding  temperature  ranges  were  308  -  720  K  and  296  -  520  K, 
respectively.  Room  temperature  measurements  were  also  made  In  argon  for  both 
Na  and  K.  Pressures  encompassed  at  each  temperature  range  from  1  to  between  5 
and  8  torr.  Above  700  K,  no  reliable  decays  could  be  measured,  because  alkali 
atoms  which  had  accumulated  on  the  walls  at  lower  temperatures  began  to  effuse 
into  the  flow,  interfering  with  the  rate  measurements. 

Using  the  data  for  Na  +  O2  +  (He)  as  an  example,  we  see  in  Fig.  2  that 
the  first-order  (logarithmic)  decays  of  signal  with  reaction  time  are  linear 
for  over  a  factor  of  100.  Plots  of  a  series  of  first-order  rate  constants 
(corrected  for  wall  and  dlffuaional  effects)  versus  [02]  are  illustrated  In 
Fig.  3.  Second-order  rate  constants  derived  from  these  are  then  plotted  vs. 
total  number  density  and  the  third-order  rate  constants  are  determined  from 
the  slopes  of  these  plots.  In  most  cases,  these  lines  Intersect  the  origin  to 
within  the  statistical  uncertainties  of  the  fit  (±10  s-1).  For  Na  +  02  + 

(N2),  a  small  positive  intercept  is  observed  at  all  temperatures.  The  cause 
for  the  Intercept  is  unclear.  An  explanation  sometimes  given  for  such 
behavior  is  that  the  wall  acts  as  a  stabilizing  third  body  and  adds  a  pressure 
Independent  component  to  k**.  However,  this  explanation  would  appear 
unlikely  in  these  experiments  given  the  near  unit  efficiency  for  wall  removal 
of  all  alkali  species.  An  alternative  explanation  is  that  insufficient 
diffusional  mixing  of  the  injected  O2  could  lead  to  underestimation  of  k11 
at  higher  pressures,  resulting  in  an  apparent  intercept  in  the  k**  vs. 
pressure  plots.  However,  calculations  of  mixing  length  indicate  that  this 
effect  should  not  be  significant  even  at  the  highest  pressure  (8  torr)  in 
these  studies. 
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The  three-body  association  rate  constants  as  a  function  of  temperature 
for  various  third  body  collision  partners  are  shown  for  sodium  in  Fig.  4  and 
for  potassium  in  Fig.  5.  The  data  can  be  described  by  an  expression  of  the 
form  k***(T)  -  AX“n.  The  results  of  nonlinear  least-square  fits  to  this 
function  are  found  in  Table  I.  For  sodium  in  helium  and  nitrogen,  and  for 
potassium  in  helium,  a  simple  T~^  representation  could  also  be  used.  However, 
the  results  for  potassium  in  nitrogen  show  a  somewhat  weaker  temperature 
dependence,  best  represented  by  x“(0.56  ±0.20)  although  this  difference  is 
not  significant  within  the  statistical  uncertainties  of  the  data. 

The  expressed  total  experimental  uncertainties.  Including  allowance  for 
systematic  errors,  can  be  estimated  as  the  square  root  of  the  sum  of  the 
squared  individual  uncertainties  due  to  a)  flow  velocity,  temperature, 
pressure,  and  concentration  calibration  factors,  5 Z,  b)  the  precision  in 
determining  the  coefficient  A  in  the  expression  for  k^^CT),  17Z,  and  c)  the 
uncertainty  in  the  correction  of  observed  first-order  decays  for  wall  removal 
and  diffusional  effects,  10Z. 

Table  I  also  contains  other  determinations  of  k2.  Husain  and 
Plane28-29  used  a  flash  photolysis  system  to  generate  alkali  atoms  from  KI  and 
Nal  and  followed  the  decay  of  K  and  Na  by  atomic  resonance  absorption.  They 
observed  no  temperature  dependence  over  their  limited  temperature  range  (724  K 
and  844  K)  for  the  Na  +  02  +  M  reaction.  A  substantially  larger  variation 
with  temperature  was  obtained  for  their  K  +  02  +  M  study  although  they 
arbitrarily  fit  their  data  to  a  T-*  dependence.  The  agreement  between  their 
values  and  the  slightly  extrapolated  temperature  dependent  rate  constants  for 
both  reactions  from  this  work  is  satisfactory,  especially  considering  the  very 
different  methods  employed.  This  comparison  is  shown  in  Figs.  4  and  5. 

Hynes  et  al.25  have  studied  the  Na  +  02  +  M  reaction  in  experiments  using 
laser  Induced  fluorescence  to  measure  Na  and  OH  profiles  under  fuel  lean 
conditions  in  H2/02/N2  flames.  A  detailed  kinetic  model  assuming  reasonable 
values  for  rate  constants  for  the  reactions  interconnecting  the  species  Na, 
Na02 ,  NaO  and  NaOH  i3  used  to  determine  the  best  fit  to  the  observations. 
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Their  optimum  value  of  k(Na  +  O2  +  M),  assuming  a  T"1  temperature  dependence. 
Is  shown  In  Fig.  4,  and  is  lower  by  about  a  factor  of  4  compared  to  the 
extrapolated  fit  from  this  work.  Part  of  this  difference  may  be  explained  by 
the  different  identities  of  M  (N2,  02,  H20  mixtures)  in  their  flame  studies. 
-Hynes  et  al.25  also  explain  the  discrepancy  between  their  result  and  the 
earlier  flame  studies, 17-19  which  obtained  a  lower  value  for  k(Na  +  02  +  M)  by 
nearly  three  orders  of  magnitude.  The  earlier  works  had  incorrectly 
attributed  the  global  disappearance  of  Na  to  reaction  with  02,  while  in  fact 
the  Na  concentration  follows  that  of  atomic  hydrogen,  whose  decay  is 
controlled  by  reaction  with  02  to  form  H02.  A.  similar  explanation  has  been 
proposed  by  Jensen.24 

The  much  earlier  diffusion  flame  studies15"1,6  correctly  identified  the 

Na  +  02  +  M  reaction  tp  be  fast  for  a  three  body  process.  The  later  and  more 

extensive  study  by  Bawn  and  Evans16  obtained  a  value  at  533  K  which  is  higher 

by  a  factor  of  3  than  this  study.  They  also  observed  a  fall-off  behavior  for 

their  second  order  rate  constant  with  pressure  above  about  10  torr  and  present 

their  results  in  terms  of  an  energy  transfer  mechanism  using  a  Lindemann-type 

expression  to  interpret  their  data.  This  fall-off  behavior  was  not  observed 
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by  Husain  and  Plane*  who  show  linear  plots  of  k2  versus  pressure  from  25  to 
150  torr.  Their  effective  bimolecular  rate  at  150  torr  and  724  K  for  the  Na  + 
02  +  (N2)  reaction  of  2.2  x  10"12  cm3  s"1  is  nearly  equal  to  the  high  pressure 
limit  inferred  by  Bawn  and  Evans  of  3  x  10" 12  cm3  s"1.  It  thus  appears  that 
the  pressure  dependence  observed  by  Bawn  and  Evans  may  have  been  an  artifact 
of  their  experimental  diffusion  flame  technique.  No  curvature  was  observed  in 
our  lower  pressure  P  <  8  torr  conditions  of  this  study  and  under  these 
conditions,  the  recombination  reaction  appears  to  be  in  its  low  pressure, 
third  order,  limit. 

In  this  limit,  the  recombination  rate  constant  can  best  be  calculated 
using  unimolecular  rate  theory,34  In  which  second-order  rate  constants  for 
unimolecular  decomposition  are  related  through  the  equilibrium  constant  to  the 
third  order  (low  pressure)  association  rate  constants.  The  unimolecular  rate 
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constants  are  obtained  using  simplified  equations  of  Troe  based  on  RRKM 

theory.35-36  The  recombination  rate  constant  can  be  expressed  as  the 

sc 

product  of  a  strong  collision  rate  constant  kfec  q  and  a  weak  collision 
deactivation  efficiency  term  gc> 
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In  effect,  8C  is  a  term  expressing  the  efficiency  by  which  (Na02)*  is 
stabilized  upon  collision.  In  the  absence  of  detailed  state-to-state  energy 
transfer  rates,  8c  is  an  adjustable  parameter  which  depends  on  the  average 
change  in  internal  energy  of  the  transition  state  complex  per  collision,  <aE>, 
so  that 


<AE> 

f£rt 


(7) 


Fg  is  a  correction  term  for  the  energy  dependence  of  the  density  of 
states  in  the  transition  complex.  Using  a  value  of  163  kJ  mole-1  for  the 
Na-02  bond  energy,25  Patrick  and  Golden37  have  calculated  k^gC  for  Na  +  02  + 

SC 

(He)  at  300  and  700  K.  The  computations  result  in  krec  (300  K)  *  7.8  x 

lO-30  cm6  s-1  for  M  ®  He.  Thus,  a  value  of  8C  *  0.18  is  required  to 

reproduce  the  experimental  result.  The  average  energy  transferred 

per  collision  is  0.8  kJ  mole-1 .  At  700  K,  kSC  *  4.2  x  10-30  cm6 

rec,o 

s-1  ,  suggesting  a  value  of  8c  *  0.09  if  <AE>  is  assumed  to  be  independent  of 
temperature.  Extrapolation  of  the  experimental  value  for  helium  to  700  K 
results  in  an  observed  8C  =*  0.14,  in  good  agreement,  considering  the 
simplified  assumptions  of  this  approach. 

This  unimolecular  approach  has  been  shown  to  be  fairly  reliable  in 
calculating  termolecular  association  rate  constants  for  a  variety  of 
species.34*36  The  values  of  8C  for  M  =  N2  and  Ar  tend  to  lie  in  the  range 
of  0.1  -  0.5.  The  calculations  for  Na  +  02  in  helium  agree  quite  well  with 
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these  numbers,  given  that  one  might  expect  helium  to  be  somewhat  less 
efficient  than  N2  in  quenching  the  excited  intermediate.  For  the  similar 
reaction  H  +  02  +  (N2),  the  calculations  also  perform  well,  giving  0c  “ 

0.08,  even  though  the  rate  constant  for  this  reaction,  S.9  x  10“32  cm6  s_1  at 
300  K,38  is  approximately  30  times  slower  than  the  corresponding  sodium 
reaction. 

In  extrapolating  our  measured  rate  constants  to  higher  pressures,  it  is 
important  to  understand  the  fall-off  behavior  as  the  reaction  mechanism  goes 
from  its  low  pressure  third  order  limit  to  its  high  pressure  second  order 
limit.  A  method  for  calculating  the  fall-off  parameter  based  on  RRKM  theory 
has  been  developed  by  Luther  and  Troe.39  Using  their  approach,  the 
bimolecular  rate  constant  in  the  transition  region  may  be  calculated  as34*38 


II 


k.  +  k0tKJ  ? 


(8) 


where  and  Icq  are  the  high  pressure  and  low  pressure  limiting  rate 
constants,  respectively.  The  broadening  parameter,  F,  is  dependent  on  the 
molecular  structure  of  the  adduct  and  typically  is  on  the  order  of  0.6  for 
small  molecules. 3 4  The  Patrick  and  Golden  calculation  for  Na02*  gives  F  - 
0.5. 37 

In  order  to  assess  the  fall-off  behavior,  the  association  rate  k„ 
may  be  estimated  from  the  dynamics  of  the  Na  +  02  encounter.  Alkali  atom 
reactions  are  often  described  in  terms  of  their  ionic  character  which  leads  to 
very  fast  reaction  rates  via  an  electron  jump  mechanism.40  Ionic  forces  in  the 
alkali-oxygen  associations,  however,  are  not  as  dominant  as  in  other  alkali 
atom  reaction  mechanisms.  Although  the  approaching  species  do  adiabatically 
transfer  from  the  incident  covalent  potential  surface  to  an  ionic  surface 
(electron  jump  mechanism),  the  crossing  distance  (rc)  is  relatively  small  as 
compared  with  alkali-halogen  systems  and  the  overall  association  rate  is 
dominated  by  the  dispersion  interaction  on  the  covalent  potential  surface. 


This  Is  Illustrated  by  comparing  the  rate  constant  based  on  the  product  of  the 
mean  collision  velocity  and  the  electron  jump  cross  section  (modified  to 
account  for  the  dispersion  and  angular  momentum  terms  in  the  covalent 
potential  function),4*  to  a  rate  constant  computed  solely  from  a  close 
collision  model,42  in  which  all  collisions  surmounting  the  angular  momentum 
barrier  of  a  Cr“6  potential  contribute  to  the  rate.  For  Na  +  O2  at  300  K,  an 
electron  jump  rate  constant  (rc  ■  2.6  A)  is  1.5  x  10"10  cm3  s”1,  while  the 
close  collision  constant  is  5.9  x  10” 10  cm3  a"1.  Thus,  it  appears  that  the 
rate  of  Intermediate  complex  formation  for  reaction  (2)  is  governed  by  long 
range  dispersion  forces,  as  found  in  most  covalent  systems.  Once  formed, 
however,  these  intermediates  are  certainly  (alkali)+02”  ion  pairs.22  In 
either  case,  the  association  rate  constant  is  expected  to  be  fast  with  a  value 
for  k.  >  10“10  cm3  a”1. 

Using  a  value  for  k«  -  2  x  10"1Q  cm3'  s-1,  Eq.  (8)  may  be  used  to 
estimate  the  ratio  of  the  observed  rate  constant,  k^g,  to  the  true  third 
order  low  pressure  limit,  kg.  At  the  maximum  density  in  our  studies,  1.6  x 
1017  cm”3,  one  obtains  kobs/^o  -  0.93  which  Indicates  that  our  experiments 
are  essentially  in  the  low  pressure  limit.  For  the  conditions  of  Husain  and 
Plane28  (150  torr,  724  K) ,  this  same  calculation  would  predict  kgb8/k0  - 
0.87.  For  the  combustion  conditions  with  P  -  760  torr  T  -  2000  K  applicable 
to  the  Hynes  et  al.  study,25  one  obtains  kgb8/k0  ■  0.89.  Thus,  the  large 
value  for  k*,  would  indicate  that  both  of  these  studies  are  also  essentially 
in  the  low  pressure  limit  with  regard  to  the  recombination  reaction. 
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FIGURE  CAPTIONS 


Figure  1 
Figure  2 


Figure  3 


Figure  4 


Figure  5 


Schematic  view  of  the  high  temperature  flow  reactor. 

Pseudo-first  order  decay  of  Na  fluorescence  signal  vs.  reaction 
distance  for  M  -  He  at  309  K  and  5.05  torr.  [o2]  ■  1.09  (a),  3.36 
(b),  7.29  (c),  14.3  (d),  26.6  (e);  units  of  1014  cm-3.  The 
estimated  relative  standard  deviations  for  these  points  range  from 
1Z  to  4Z. 

Plots  of  first-order  Na  decay  rates  vs.  [02]  for  M  ■  He  and  T  ■ 
309-335  K.  Total  pressure  (torr)  •  1.05  (o),  1.51  (e),  2.52 
(+),  3.04  (a),  5.05  (7).  Each  point  has  been  corrected  for  the 
effects  of  diffusion  and  wall  loss.  Error  bars  correspond  to  ±2 a. 
Plot  of  In  k(T)  vs.  In  T  for  the  reaction  Na  +  02  +  M,  with  M  -  N2 
(circles).  He  (triangles),  and  Ar  (square)  for  this  work.  The 
solid  lines  are  least-square  fits  to  these  data.  The  +  and  x 
symbols  represent  data  for  M  ■  N2  and  He,  respectively  from 
reference  28,  and  the  dashed  line  (M  -  flame  gases)  Is  from 
reference  25. 

Plot  of  In  k(T)  vs.  In  T  for  the  reaction  K  +  02  +  M,  with  M  ■  N2 
(circles).  He  (triangles),  and  Ar  (square).  The  solid  lines  are 
least-square  fits  to  these  data.  The  4-  and  x  symbols  represent 
data  for  M  ■  N2  and  He,  respectively  from  reference  29. 
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APPENDIX  B 


THE  GAS  PHASE  REACTION  RATE  OF  SODIUM  HYDROXIDE  WITH  HYDROCHLORIC  ACID 
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The  reactions  of  metallic  species  introduced  into  the  atmosphere  by  meteor  ablation  may  play  a  significant  role  in  stratospheric 
chemistry.  In  particular,  it  has  been  suggested  that  the  reaction  of  NaOH  with  HC1  might  affect  the  concentration  of  odd 
chlorine,  thus  having  an  impact  on  the  ozone  balance.  This  paper  describes  the  first  measurement  of  this  reaction  rate  consiam 
At  308  K,  we  find  that  k  *  (2.8  ±  0.9)  x  10"10  cmJ  molecule'1  s'1  As  a  result  of  the  methods  developed  to  perform  this 
measurement,  we  have  also  determined  estimates  of  the  following  room  temperature  rate  constants  in  units  of  cm’  molecule  ’ 
s  ':  A(NaO  +  HCI  —  NaCl  +  OH)  =*  2.8  x  10"'°,  K(NaOH  +  H  —  Na  +  H20)  >  4  x  t<r,a.  *<NaCI  +  H  -  Na  +  HCI) 
a  5  x  10^'**”,  and  Ic(Na  +  H202)  »  6.9  x  I0'n,  where  approximately  0.6  of  the  reactions  produce  NaOH  4  OH.  with 
the  remainder  forming  NaO  +  H20. 


Introduction 

Metallic  elements  volatilized  during  meteor  entry  into  the 
Earth's  upper  atmosphere  play  a  significant  role  in  the  structure 
of  the  D  and  H  regions  of  the  ionosphere,1'2  and,  at  least  in  the 


(1)  E.  Murad,  J  Geophys  Res  .  *3.  5525  (1978) 


case  of  sodium,  the  visible  day  and  ntghtgiow  emissions  from  the 
mesosphere  and  lower  thermosphere.2'’  Recently,  it  was  suggested 


(2)  T  L  Brown,  Chem.  Rev  ,  73,  64S  (1973) 

(3)  1  E  Blamont  and  T  M  Donahue.  J  Geophys  Res .  69  4093  ( l»S4) 

(4)  D  M  Hunten,  Space  Set  Rev  ,  6.  493  1 1967) 
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Figure  1.  Schematic  diagram  of  atmospheric  sodium  chemistry 


that  sodium,  and  other  meteor  metals  may  be  important  in 
stratospheric  chemistry  by  affecting  ozone  reduction  by  the 
catalytic  chlorine  cycle.4-* 

The  influx  of  meteor  metals  into  the  upper  atmosphere  has  been 
estimated9  to  be  3.5  x  104  kg  yr'1,  with  a  sodium  abundance  of 
2%  leading  to  a  calculated  sodium  flux  of  1.2  x  (04  atoms  cm'3 
s'1.  Other  estimates  of  sodium  flux  run  as  high  as  2  X  104  cm'2 
s'1.10  The  flux  of  other  metallic  species  such  as  Mg,  Ca,  Al,  Si, 
and  Fe  will  be  as  much  as  10  times  higher  and  speculative  concerns 
about  their  influence  on  upper  atmospheric  homogeneous  and 
heterogeneous  chemistry  have  been  published  1'2-,'*1° 

Unfortunately,  all  attempts  to  model  the  role  of  volatilized 
meteor  metals  (particularly  sodium)  in  the  mesosphere  and 
stratosphere3"* 7  9-15  have  suffered  from  an  almost  total  lack  of 
measured  rate  constants.  All  such  models  start  with  the  oxidation 
of  sodium  or  other  metallic  species  in  reaction  with  atmospheric 
O,  Oj,  or  O,.  However,  the  only  measured  chemical  rate  constants 
currently  available  for  any  alkali  meteor  metal  oxidation  reactions 
are  those  for  the  three-body  recombination  reactions  of  alkali 
atoms  (Na,  K)  with  CV4-1* 

Na  or  K  +  02  +  M  -  Na02  or  KQ2  +  M  ( 1 ) 


(5)  (a)  V  W  J  H  Kirchhoff,  B  R  Clemeatu,  and  D  M  Simomch.  ] 
Geophys  Res  .  M.  1323  (1979),  (b)  D  R  Bales  and  P  C  Ojha.  .Valuer 
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(8)  E  Murad.  W  Swider.  and  S  W  Benson.  Mature  (London).  289  273 
(1981) 
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i  12)  V  W  J  H  Kirchhoff,  B  R  Clemesha.  and  D  M  Simomch.  J 
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In  fact,  until  2  years  ago,  even  this  reaction  was  thought  to  be 
— 1000  times  slower19-20  than  the  value  found  in  the  more  recent 
measurements. 

From  the  recent  modeling  work7'*-10-13-21'22  and  comparison  of 
these  models  with  atmospheric  measurements,12-14-13  neutral  so¬ 
dium  is  believed  to  be  transformed  via  a  series  of  chemical  re¬ 
actions  involving  NaO  and  Na02  intermediates  to  NaOH.  A 
survey  of  the  relevant  literature4-*-1013  provides  a  fairly  complete 
list  of  possible  neutral  sodium  reactions  of  importance,  and  a 
simplified  schematic  diagram  of  these  reactions  is  shown  in  Figure 
I. 

Neutral  sodium  appears  in  a  range  between  1 10-  and  70-km 
altitude.  As  shown  in  Figure  1 ,  the  main  removal  mechanisms 
are  by  reaction  with  02  or  0},  the  latter  used  by  Chapman23  to 
explain  the  Na  nightglow.  We  have  recently  completed  tem¬ 
perature-dependent  measurements  of  the  reaction  rates  with  O, 
in  our  laboratory  which  show  that  this  reaction  proceeds  with  a 
fast  three-body  rate  constant  of  1.9  x  10-30  cm6  s-r  with  N2  as 
the  third  body  at  300  K."  This  leads  us  to  believe  that  most  of 
the  Na  in  the  lower  mesosphere  is  converted  to  Na02,  since  the 
three-body  recombination  rate  with  02  exceeds  even  a  gas  kinetic 
two-body  rate  with  O3  below  80  km.  Since  none  of  the  remaining 
rates  have  ever  been  measured,  the  rest  of  the  mechanism  (as 
proposed  by  various  models)  is  speculative.  However,  the  dominant 
sodium  species  which  leave  the  mesosphere  and  enter  the  stra¬ 
tosphere  are  probably  NaOH  and/or  Na02. 

The  fate  of  NaOH  is  very  uncertain.  Indeed,  only  in  the  past 
few  years  has  NaOH  been  recognized  as  a  major  reaction  product. 
Ferguson4  suggested  that  it  forms  NaOH  cluster  ions  of  the  form 
H*(NaOH),,(H20)„,  which  either  may  be  rained  out  from  the 
troposphere  or  removed  in  the  stratospheric  aerosol  layer.  In  a 
recent  paper  by  Murad  et  al..*  it  was  proposed  that  the  reactions 
of  metal  hydroxides  (and  superoxides)  with  chlorine  compounds 
between  40  and  70  km  may  have  an  impact  on  the  depletion  of 
stratospheric  ozone.  In  the  case  of  sodium,  the  exothermic  bi- 
molecular  reactions 

NaOH  +  HCI  -  NaCl  +  H20 
iff,  -  -133  ±  13  kJ  mol-1 

NaOH  +  CIO  —  NaCl  +  H02 
UP",  -  -63  ±  15  U  mol-1  (3) 

NaOH  +  Cl  —  NaCl  +  OH 

«  -67  ±  13  kJ  mol'1  (4) 

might  be  expected  to  proceed  rapidly  and  act  as  a  sink  for  Cl. 
given  that  NaCl  can  readily  polymerize  and  condense  via  het¬ 
erogeneous  nucleation.'  Murad  et  al  calculated  thal  if  Jfc2  were 
~  10  "  cm3  molecule"1  s'1,  then  the  reaction  of  NaOH  with  HCI 
would  be  comparable  to  the  major  recognized  Cl  regeneration 
mechanism 

OH  +  HCI  —  Cl  +  H20  (5) 

While  previously  published  studies  have  viewed  NaCl  as  a  potential 
sink  for  stratospheric  chlorine,7'*  more  recent  analyses  by  Row- 
'  land24  indicate  that  photolysis  of  NaCl  may  in  fact  release  free 
Cl  Given  the  potentially  large  J  values  (photolysis  rates)  for  this 
process,24  reactions  2-4  could  effectively  supplement  reaction  5 
as  a  release  mechanism  for  Cl  from  the  inactive  HCI  stratospheric 
reservoir  and  thereby  determine  the  extent  to  which  ozone  might 
be  depleted  by  chlorine  compounds  in  the  stratosphere  However. 

R  Carabetta  and  W  E  Kaikan.  J  Phix  Chem  .  72.  ^483  ( J 9*8) 
(20)  M  J  Met  wan  and  l.  F  Phillip*.  Trans  Faradav  S*  jc  .  62,  7 1  7 
(1966) 
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(22)  D  M  Simomch.  B  R  Clemesha.  and  V  W  J  H  Kirchhoff,  J 
OtopMvi  Rex  *4  I <43  (  I  979) 
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Figure  2.  Schematic  view  of  the  flow  tube 


to  understand  fully  the  role  of  alkali  species  in  the  atmosphere, 
one  must  also  consider  the  effects  of  Na02  and  NaOH  photo¬ 
dissociation  on  these  processes. I3J3 

The  purpose  of  this  study  is  to  provide  a  direct  experimental 
rate  measurement  of  the  reaction  NaOH  +  HC1  —  NaCl  +  H2O, 
which  serves  as  a  starting  point  for  understanding  the  stratospheric 
role  of  sodium  and  other  meteor  metals.  Only  by  obtaining 
accurate,  directly  measured  rate  constants  can  we  hope  to  un¬ 
derstand  atmospheric  metal  chemistry,  and,  in  particular,  how 
this  chemistry  affects  the  ozone  balance. 

Experimental  Section 

Very  little  gas-phase  research  has  been  done  on  alkali  hydroxides 
because  they  are  difficult  to  vaporize,24  are  extremely  corrosive,21 
and  readily  dimerize  in  the  gas  phase. 21  In  addition,  there  are 
no  established  detection  techniques  sensitive  enough  to  allow 
kinetic  analyses  of  alkali  hydroxide  reactions.  With  this  in  mind 
we  have  taken  an  indirect  approach  for  producing  and  detecting 
NaOH,  i.e.,  chemical  production  of  NaOH  by  the  reaction  of 
atomic  sodium  with  hydrogen  peroxide,  and  detection  of  NaOH 
by  chemical  conversion  back  to  atomic  sodium,  which  is  then 
observed  by  laser-induced  fluorescence. 

The  measurements  were  performed  in  the  Aerodyne  high- 
temperature  fast-flow  reactor,  which  is  fully  described  in  ref  29 
and  whose  relevant  features  are  illustrated  in  Figure  2.  Briefly, 
a  7.26-cm-diamcter,  1 20-cm-long  alumina  tube  is  used,  fitted  with 
Kanthal  heater  elements  which  can  radiatively  heat  the  tube  over 
the  temperature  range  294-1500  K.  Four  perpendicular  alumina 
side  arms  at  the  tube  exit  permit  detection  of  flow  species  by  a 
variety  of  optical  techniques,  including  laser-  or  resonance 
lamp-excited  fluorescence,  infrared  absorption  (either  broad  band 
or  high  resolution),  and  observation  of  chemiluminescence.  Gas 
temperatures  are  obtained  with  chromei-aJumel  or  shielded  W-5% 
Re — W-26%  Re  thermocouples.30  Extensive  calibrations  indicate 
the  temperature  can  be  measured  with  an  accuracy  of  ±10  K  over 
the  entire  operating  range  of  the  reactor  for  flow  Reynolds 
numbers  below  500;  for  the  current  experiments  they  are  typically 
below  50.  The  gas  is  pumped  by  a  Kinney  KMBD  1602  me¬ 
chanical  pump  and  Roots  blower  with  an  effective  pumping  speed 
of  450  L  s'1.  The  helium  carrier  gas  is  added  at  the  entrance  of 
the  flow  tube  through  mullite  multichannel  arrays  which  lami- 
narize  the  flow.  These  are  ~2.5-cm  upstream  of  the  outlet  of 


(25)  f  S.  Rowland  and  Y  Makide.  Geophys  Res  Lett.,  9,  473  (1982). 
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(29)  M.  E.  Gerah,  J  A.  Silver.  M.  S.  Zahniser.  C  E  Kolb.  R  G  Brown, 
C  M.  Gozewaki,  S  Kallelis.  and  J  C  Wormhoudt,  Rev  Set.  Instrum.,  52. 
1213  (1981) 
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the  reactant  inlet  tube  and  78  cm  from  the  detection  region.  This 
allows  sufficient  distance  ( 1 8  cm)  for  the  helium  flow  to  mix  with 
the  reactants  and  develop  a  parabolic  velocity  profile  before 
reaching  the  reaction  zone.  Gas  volumetric  flow  rates  are  de¬ 
termined  with  calibrated  thermal  conductivity  type  mass  flow 
meters.  Flow  speeds  can  be  varied  from  4  to  100  m  s'1  A 
calibrated  MKS  Baratron  Model  310-BSH10  capacitance  ma¬ 
nometer  (0.8%  accuracy)  is  used  to  measure  pressure. 

Alkali  atoms  are  generated  by  heating  the  sample  in  a  2.5- 
cm-diameter  cylindrical  monel  oven  to  a  temperature  commen¬ 
surate  with  attaining  a  vapor  pressure  of  that  species  of  1 O'*  to 
10~*  torr  within  the  oven.  The  oven  is  silver  plated  to  resist  alkali 
corrosion.27  The  vapor  is  entrained  in  a  flow  of  inert  carrier  gas 
and  introduced  into  the  flow  tube  through  a  10-cm  section  of 
19-mm  o.d.  silver  tubing.  The  sodium  vapor  is  further  diluted 
by  the  carrier  gas  in  the  main  flow  tube  so  that  the  sodium 
concentration  within  the  reaction  zone  is  always  less  than  1010 
cm'3.  For  wall  removal  rate  measurements  the  entire  oven  as¬ 
sembly  may  be  placed  downstream  of  the  mullite  arrays  as  a 
movable  source.  Since  the  oven  is  heated,  it  warms  the  main 
carrier  gas  flow  slightly  Axial  temperature  surveys  in  the  reaction 
region  show  that  the  final  flow  temperature  profile  is  uniform  at 
a  value  of  308  K. 

Sodium  hydroxide  is  produced  via  the  reaction 
Na  +  H202  —  NaOH  +  OH 
Mf°r  *=  -130  ±  13  kJ  mol'1 

The  H202  is  added  upstream  of  the  mullite  arrays  with  the  main 
carrier  gas.  An  all  glass  and  teflon  inlet  system  is  used  to  prevent 
decomposition  of  H202  between  its  reservoir  and  the  flow  tube 
Flow  rates  of  H202  are  determined  by  entraining  the  vapor  in  a 
measured  helium  flow  at  a  known  total  reservoir  pressure,  where 

H202  flow  =»  (7>H,o,/J>!oui)(He  flow  through  system)  (7) 

The  partial  pressure  of  H202  at  25  °C  is  2.0  torr. 

Detection  of  NaOH  is  accomplished  by  converting  it  back  to 
atomic  sodium  in  the  detection  region,  where  the  sodium  is 
measured  by  laser-induced  fluorescence  (L1F).  This  conversion 
is  accomplished  by  injecting  an  excess  of  atomic  hydrogen  into 
the  flow  2-cm  upstream  of  the  LIF  detector. 

NaOH  +  H  —  Na  +  H20 

(  $  \ 

AH°r  »  -154  ±13  kJ  mol'1 

Na  +  hv  -»  Na*  -*  Na  +  hv  LIF  detection  (9) 

The  hydrogen  atoms,  produced  by  microwave  discharge  of  pure 
molecular  hydrogen,  are  introduced  through  a  teflon-lined  6- 
mm-diameter  tube  With  the  mean  flow  velocity  of  1000  cm  s 
the  reaction  time  of  hydrogen  in  the  detector  (rd)  is  2  ms  The 
atomic  hydrogen  concentration  is  estimated  to  be  2  x  10“  cm  ’, 
by  measuring  the  H2  flow  rate  and  assuming  that  10%  of  the  H; 
passing  through  the  discharge  dissociates,31  with  no  recombina¬ 
tion32  or  loss  on  the  inlet  tube  walls  31 
A  Molectron  DL14  nitrogen  pumped  dye  laser  is  used  foi 
laser-induced  fluorescence  detection  of  Na.  The  laser-induced 
fluorescence  data  acquisition  system  and  manipulation  of  data 
have  been  detailed  elsewhere.33  However,  it  should  be  noted  that 
fluorescence  is  usually  averaged  over  100  laser  pulses,  accounting 
for  nonfluorescent  background  signals  and  for  pulse-to-pulse 
fluctuations  in  laser  intensity.  The  combined  signal-to-noise  ratio 
for  these  measurements  generally  exceeds  25 

Purities  of  the  chemicals  used  in  these  experiments  are  as 
follows,  sodium  metal.  99.95%  (Alfa);  helium.  99  995% 
(Northeast  Cryogenics),  hydrogen,  99  995%  (Air  Products),  ni¬ 
trogen.  99.998%  (Northeast  Cryogenics);  and  hydrogen  chloride. 
99  99%  (Northeast  Cryogenics)  Hydrogen  peroxide,  obtained 


(31)  W  E  Jones,  S  D  MacKnight,  and  L  Teng,  Chem  Rev.  73.407 
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TABLE  1:  Typical  Experimental  Condition 


temp,  K 

flow  velocity,  m  s'1 
press.,  tori 
[Na],  cm’1  (initial) 
[HA),  cm’1 
(HCll,  cm"’ 

(HI.  cm'1 


308 
10.0 
2.0 
SIO10 
1.3  X  10” 
(1-4)  X  10” 
2  X  1014 


as  a  90%  (by  weight)  solution  from  FMC  Corp.,  is  purified  by 
vacuum  distillation.  Titration  with  potassium  permanganate 
indicates  a  resulting  purity  on  the  order  of  99%.  This  corresponds 
to  a  purity  in  the  vapor  phase  of  ~93%,  with  the  balance  being 
water  vapor. 

Rate  measurements  are  made  with  one  reactant  in  excess  of 
the  other  (detected)  reactant,  thus  ensuring  pseudo-first-order 
kinetic  conditions.  Reaction  times  are  varied  by  changing  the 
injector  position  at  fixed  total  flow  velocity  and  pressure.  Cor¬ 
rections  for  both  axial  and  radial  diffusion  and  wall  removal  are 
made  with  the  procedure  outlined  by  Brown.54  This  method  is 
based  on  a  numerical  solution  of  the  equations  describing  diffusion 
and  reaction  in  a  flow  tube.  It  assumes  that  PoiseuiUe  flow  exists 
and  provides  ranges  for  k,  k,  (wall  removal),  and  D  (diffusion 
coefficient)  for  which  the  solutions  have  been  shown  to  be  valid. 
Wall  removal  rates  of  the  measured  species  are  determined  in 
separate  experiments  by  varying  the  oven  position.  The  observed 
wall  removal  rate,  as  well  as  the  reaction  rate,  is  corrected  for 
diffusion  effects.  Diffusion  constants  for  alkali  atoms  and  hy¬ 
droxides  were  obtained  from  wall  removal  measurements  in  the 
instance  where  the  observed  disappearance  of  the  species  is  dif¬ 
fusion  limited.”  At  308  K,  />(Na-He)  *  0.50  atm  cm2  s"1  and 
D(NaOH-He)  *  0.47  atm  cm2  s'1. 

An  important  aspect  of  the  data  analysis  is  the  ability  to  ensure 
that  secondary  reactions  have  no  effect  on  the  concentration  of 
the  species  being  monitored.  This  could  dramatically  affect  the 
accuracy  of  the  rate  measurements.  The  effects  of  secondary 
reactions  are  determined  by  modeling  the  reactions  occurring  in 
the  flow  tube.  This  is  done  with  the  Aerodyne  package  code,5* 
a  kinetic  modeling  program  which  numerically  integrates  the 
differential  rate  equations  for  a  specified  set  of  reactions. 
Backward  reaction  rates  determined  from  the  JANAF  thermo- 
chemical  tables  and  the  forward  rates  are  included,  to  ensure  that 
accuracy  is  maintained. 


Results 

Na  +  //A  Reaction.  The  NaOH  formation  rate  was  measured 
by  directly  observing  the  disappearance  of  sodium,  with  H202  in 
known  excess.  The  reaction  has  two  exothermic  product  channels 
Na  +  H,0,  —  NaOH  +  OH 

(6) 


A//°r«-130±  13  kJ  mol"1 
Na  +  HjO,  —  NaO  +  HjO 
\H°,  -  -130  ±  42  kJ  mol'1 


(«') 


so  what  is  measured  by  monitoring  Na  disappearance  is  the  total 
reaction  rate.  However,  one  can  add  excess  CO (~2  x  10”  cm"5), 
which  rapidly  and  quantitatively  converts  NaO  to  Na57'5* 

NaO  +  CO  -»  Na(2S)  +  CO, 

Atf°,  -  -259  ±  42  kJ  mol'1  (10) 


NaO  +  CO  —  Na(aP)  +  CO, 
*  -56  ±  42  kJ  mol"1 


(un 
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Flg«e  3.  Pseudo-first-order  decays  for  Na  +  H,0,  reaction.  The  solid 
line  is  the  total  reaction  rate;  the  dashed  line  is  the  rate  for  branching 
to  NaOH  only  (CO  added). 

Although  the  value  for  kw  has  not  been  measured,  indirect  es¬ 
timates  from  ref  37  imply  kl0  >  10~n  cm3  molecule-1  s'1.  The 
large  amount  of  CO  added  is  required  to  ensure  that  all  of  the 
NaO  cycles  back  to  Na  on  a  time  scale  much  shorter  than  the 
time  required  for  reaction  6  to  occur. 

A  rate  measurement  with  CO  present  results  in  the  production 
rate  for  only  the  NaOH  branch.  The  results  of  these  measure¬ 
ments  are  shown  in  Figure  3.  Experimental  conditions  are  given 
in  Table  1.  The  rate  constant  for  both  channels  is  kM.  *  (6.9 
±  3.0)  X  10-"  cm3  molecule-1  s-1,  and  the  fraction  in  the  NaOH 
product  channel  is  0.61  ±  0.10.  The  major  uncertainty  in  the  rate 
constants  kt  and  k^  is  the  H,0,  concentration,  which  is  not  directly 
measured  but  is  obtained  as  indicated  earlier. 

An  independent  determination  of  the  product  branching  ratio 
for  reaction  6  is  obtained  by  observing  the  hydroxyl  radical  via 
laser-induced  fluorescence  at  308.6  nm.  When  Na  reacts  with 
HjO,,  only  the  NaOH  product  channel  produces  OH.  However, 
when  excess  HCI  is  added  to  the  NaO  and  NaOH  product 
mixture,  additional  OH  is  formed  by  the  reaction 

NaO  +  HCI  —  NaCl  +  OH 

A//°r  *  -134  ±  42  kJ  mol-1  (U) 

Thus  a  measurement  of  the  amounts  of  OH  produced  before  and 
after  addition  of  HCI  provides  the  relative  amounts  of  NaOH  and 
NaO  originally  formed  from  Na  +  H,0,.  The  observed  NaOH 
product  fraction  by  this  method  is  0.60  ±  0.10,  in  excellent 
agreement  with  the  value  obtained  from  the  rate  measurements. 

NaOH  +  HCI  Reaction.  In  the  presence  of  excess  CO.  the 
reaction  of  sodium  with  hydrogen  peroxide  produces  only  NaOH 
With  this  reaction  used  as  a  source  for  NaOH,  a  series  of  rate 
measurements  were  made  for  the  reaction  of  NaOH  with  HCI. 
First-order  decays  were  linear  for  more  than  a  factor  of  10  in 
fluorescence  signal,  with  the  HC!  concentrations  ranging  from 
1  X  1012  to  4  X  10”  cm-1.  The  results  of  these  experiments  are 
shown  in  Figure  4,  and  the  rate  constant  for  this  reaction  is  k2 
*  (2.8  ±  0.9)  x  10-l°  cm3  molecule-1  s-1.  The  uncertainty  ex¬ 
pressed  includes  estimated  precision  errors  at  a  95%  confidence 
level,  as  well  as  estimated  errors  in  accuracy 
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HCI  CONCENTRATION  (1013  molKulM  cm'3) 

Figure  4.  First-order  reaction  rates  vs.  (HCI)  for  the  reactions  of  NaOH 
(NaO)  +  HCI:  (O)  without  CO  added;  (V)  with  CO  added. 

Separate  measurements  were  also  made  without  adding  CO. 
In  this  case  NaO,  in  addition  to  NaOH,  is  present  and  reacts  with 
HCI  (reaction  1 1 ),  complicating  the  analysis.  Furthermore,  the 
detection  scheme,  based  on  conversion  to  Na  by  addition  of  atomic 
hydrogen,  would  not  appear  to  distinguish  between  NaO  and 
NaOH,  since  Na  atoms  are  also  produced  in  the  reaction 

NaO  +  H  —  Na  +  OH 
A//°,  ■  -155  ±  42  kJ  mol~‘  '  (12) 

Thus,  without  CO,  the  decrease  in  signal  upon  addition  of  HCI 
is  due  to  both  reactions  2  and  1 1 ,  weighted  by  the  branching  ratio 
from  reaction  6,6'.  The  results,  however,  show  the  same  decay 
with  and  without  added  CO  (see  Figure  4),  implying  that  both 
NaO  and  NaOH  react  with  HCI  at  approximately  the  same  rate. 
This  observation  is  reasonable  if  one  considers  the  hydroxyl  group 
on  NaOH  to  act  as  a  quasi-atom  in  respect  to  its  chemical  behavior 
given  identical  exothermicities  for  both  reactions. 

Separate  OH  measurements  confirm  that  atomic  hydrogen 
reacts  with  NaO  as  well  as  with  pure  NaOH.  Although  we  could 
not  measure  the  NaO  +  HQ  rate  constant  directly,  the  observation 
that  ( 1 )  OH  is  formed  in  the  NaO  +  H  reaction  (proving  the 
existence  of  NaO),  (2)  OH  is  also  formed  upon  addition  of  HCI 
to  NaO,  and  (3)  the  decay  rates  with  added  HCI  are  identical 
both  with  or  without  CO,  imply  that  NaO  reacts  with  HCI  at 
approximately  the  same  rate  as  NaOH  with  HCI. 

Detector  Corrections  and  Modeling.  Although  the  plots  of  In 
(signal)  vs.  reaction  time  are  linear  over  the  first  order  of  mag¬ 
nitude  decrease  in  signal,  at  longer  times  they  flatten  out  at  a  value 
typically  a  few  percent  of  the  initial  (zero  reaction  time)  signal 
(Figure  S).  This  effect  can  be  attributed  to  an  additional  com¬ 
ponent  to  the  Na  signal  from  the  reaction 

NaCI  +  H  —  Na  +  HCI  Ml9,  -  -21  ±  2  kJ  mol"1  (13) 

For  the  case  with  CO  present,  sodium  formed  in  the  detection 
region  by  addition  of  atomic  hydrogen  has  two  sources.  NaOH 
and  NaCI,  the  amount  of  each  depending  on  the  extent  to  which 
the  NaOH  +  HCI  -*■  NaCI  +  HjO  reaction  has  gone  to  com¬ 
pletion.  The  rate  equation  for  the  formation  of  sodium  in  the 
detector  is 

dINal/df  =*  /ts[NaOH]d(H)  +  lc,,[NaCI]d|H]  (14) 


Flgwt  S.  Typical  decay  of  fluorescence  signal  vs.  reaction  time  f-» 
NaOH  4-  HCI.  showing  contribution  to  total  detected  sodium  ( — )  from 
NaOH  ( - )  and  NaCI  (••■•).  Circles  are  experimental  data. 

where  the  d  subscript  indicates  the  concentration  at  the  entrance 
to  the  detector  zone.  Integrating  over  the  detector  reaction  time 
Vt),  assuming  excess  hydrogen,  we  obtain 

(Na]  -  [NaOH]d(l  -  e'WHlu)  +  [NaCl]d(l  -  e'‘'>lH"*)  (15) 

If  no  HCI  is  added,  the  observed  LIF  signal  arises  solely  from 
NaOH.  We  have  observed  that  this  signal  (after  correcting  for 
differences  in  diffusion  and  wall  removal  between  Na  and  NaOH) 
is  identical  with  that  arising  from  only  sodium  (i.e.,  without  H202 
added).  This  means  that,  within  the  available  detector  reaction 
time,  all  of  the  NaOH  is  .converted  to  sodium,  with  an  estimated 
error  of  20%.  This  sets  a  lower  limit  to  Given  /d  =  2 

ms  and  [H]  t*  2  x  1014  cm'3,  this  requires  that  lt>4x  10' 12 
cm5  molecule'1  s'1.  The  accuracy  of  tt  is  ±0.5  ms  If  less  than 
10%  of  the  H2  is  dissociated  (it  is  unlikely  to  be  higher  under  our 
operating  conditions),  then  is  faster  than  the  stated  limit. 

This  value  for  k,  is  in  reasonable  agreement  with  a  value  of 
1  X  10"IJT,/J  (1.7  x  10'"  cm5  molecule'1  s'1  at  300  K)  which 
is  used  by  Hynes  et  al.1T  in  fitting  their  flame  data  over  a  wide 
range  of  conditions,  and  also  with  2  x  iO'"<r"*00/r(l,4  x  !0'12 
at  300  K),  which  is  used  in  various  atmospheric  models. ial  3  The 
former  value  is  claimed  to  be  accurate  to  within  a  factor  of  two 
at  flame  temperatures,  and  the  latter  is  only  an  estimate.  We 
are  planning  to  perform  direct  measurements  of  this  reaction  in 
the  near  future. 

Now  if  we  measure  the  reaction  rate  of  NaOH  +  HCI,  at  short 
reaction  times  the  decay  in  signal  is  logarithmic  because  [NaOH]d 
>  (NaCl]d,  and  only  the  first  term  in  eq  IS  is  significant.  At  long 
reaction  times,  [NaCI)d  >  [NaOH]d  and  the  second  term  becomes 
dominant  (Figure  5).  If  kl}  *  It,,  we  would  never  observe  a  decay 
in  [Na]  since  both  terms  in  eq  15  would  always  have  the  same 
sum.  However,  the  fact  that  an  initial  decay  is  observed  with  a 
later  leveling  of  signal  shows  that  k  „  <<  k,  and,  fruit  the  relative 
value  of  the  signal  where  it  levels  off,  we  can  estimate  kn. 

To  determine  It,,,  we  modeled  this  reaction  system  using  the 
conditions  in  Table  11  and  rate  set  in  Table  III  This  rate  set 
is  more  extensive  than  the  above  discussion  implies  because  there 
could  be  secondary  reactions  between  the  excess  H203  and  the 
sodium  formed  in  the  detector,  as  well  as  other  minor  OH,  H, 
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TABLE  II:  laMUl  Conditions  Used  in  Chemical  Model 


temp.  K 

308 

flow  velocity,  m  s'1 

10.0 

press.,  torr 

2.0 

[Na],  cm'3  (initial) 

10'° 

[HjOj],  cm'3 

1.3  X  I0'3 

[HCI].  cm'3 

1.93  X  1011 

[HU  cm'3 

2  x  I0,3-t  x  I0'3 

lHjla.  cm'3 

10l3-10“ 

[HjO|.  cm"3 

7  X  10" 

W  ms 

0-25 

is,  ms 

0. 5-4.0 

TABLE  III:  Reaction  Rate  Set  Used  for  Modeling  Na  Chemistry* 

*,  cm3 

reaction 

molecule'1  s"1 

ref 

Na  +  H,Oj  -  NaOH  +  OH 

6  9  (-11) 

b 

NaOH  +  HCI  -  NaCl  +  H,0 

2.8  (-10) 

b 

NaOH  +  H  —  Na  +  H,0 

1.0  (-IO)-I.O  (-14) 

b 

NaCl  +  H  —  Na  +  HCI 

1.0  (-IO)-I.O  (-14) 

b 

Na  +  HOj  —  NaOH  +  O 

1.0  (-10) 

eat 

2HO,  —  HjOj  +  O, 

1.6  (-12) 

39 

O  +  HjO,  —  OH  +  HO, 

1.0  (-1  l)e'2300^r 

40 

OH  +  HjO,  —  HjO  +  HO, 

2.96  (-12)e"l44/r 

41 

H  +  HjO,  —  H,  +  HO, 

2.13  (-I2)e'l400/r 

42 

H  +  H,Oj  —  OH  +  HjO 

2.76  (— 1 2)w',*00/r 

42 

H,  +  OH  —  H  +  HjO 

7.7  (-12)«'J100'r 

40 

Hj  +  O  —  H  +  OH 

1.6  (-1 1  )e~t,n/T 

40 

OH  +  HO,  —  H,0  +  O, 

8.0  (-11) 

40 

H  +  HO,  —  20H 

3.2  (-11) 

40 

H  +  HO,  -»  H,  +  O, 

1.4  (-11) 

40 

H  +  HO,  —  H,0  +  O 

9.4  (-13) 

40 

O  +  HO,  —  Oj  +  OH 

8.0  (-1  l)r~yn,r 

40 

20H  —  O  +  HjO 

18  (-12) 

40 

HCI  +  OH  —  HjO  +  Cl 

6.6  (-13) 

40 

*  No  three-body  rates  were  used  because  the  system  is  at 

low  pres- 

jure.  ‘This  work. 


etc.  reactions.  This  model  also  includes  the  undissociated  hy¬ 
drogen.  and  the  water  impurity  (~7%)  in  the  H202.  We  sys¬ 
tematically  varied  kt,  kn,  [H],  [HjOJ.  and  rd.  The  results  show 
that  all  secondary  reactions  have  little  effect  (<5%)  on  the  cal¬ 
culated  sodium  densities  and  that  the  variations  of  Na  with  [HjOj] 
and  [H}]  are  small,  in  reasonable  agreement  with  additional 
experiments.  The  value  for  fc|3  which  best  fits  the  experimental 
data  is  (8  ±  3)/[H],  which  leads  to  value  of  k,3  in  the  range  1 
x  10'1*  to  5  x  10'13  cm3  molecule"1  s'1,  with  a  best  estimate  of 
5  x  10"'*  cm3  molecule'1  s'1.  In  light  of  the  small  exothcrmicity 
for  this  reaction,  this  value  is  not  unreasonable.  If  the  rate  constant 
is  expressed  in  the  Arrhenius  form,  k  »  Ae'tlKT,  this  would 
correspond  to  only  a  21  kj  mol"1  barrier,  even  with  a  gas  kinetic 
preexponential  term. 

Discussion 

The  reaction  rate  constant  of  HC1  with  NaOH  is  found  to  be 
in  its  gas  kinetic  limit.  Given  the  large  exothermicity,  absence 
of  obvious  steric  effects,  and  the  fact  that  reactant  and  product 
states  correlate  on  a  singlet  potential  energy  surface,  one  might 
not  consider  this  surprising.  On  the  other  hand,  the  strong  covalent 
HCI  bond  must  be  broken  during  the  reaction  with  an  efficiency 
approaching  unity.  One  can  calculate  a  rate  constant  for  this 
reaction  using  simple  close  collision  theory  43  This  procedure 
assumes  that  the  attractive  portion  of  the  intermolecular  potential 
varies  as  O'*.  where  C  depends  on  the  dipole  moments  and  static 
polarizabilities  of  the  reactants.44  and  r  is  the  intermolecular 
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distance.  The  total  effective  potential  energy  curve  at  large  r  is 
this  attractive  portion  plus  a  centrifugal  term  which  accounts  for 
orbital  angular  momentum.  The  rate  of  close  collisions  is  identified 
as  applying  to  collisions  with  total  energy  sufficient  to  surmount 
the  orbital  angular  momentum  barrier  in  this  effective  potential 
energy  curve.  If  it  is  assumed  that  all  close  collisions  react,  and 
that  there  is  no  activation  energy,  the  rate  constant  calculated 
for  NaOH  +  HCI  is  1.8  X  10"'°  cm3  molecule'1  s'*,  in  good 
agreement  with  the  measured  value. 

Few  gas-phase  reactions  of  alkali  molecules  have  been  studied, 
none  involving  alkali  hydroxides.  The  alkali-hydroxide  bond  is 
ionic  in  character,  so  that  (OH)"  is  a  closed-shell  species,  not 
dissimilar  to  its  isoelectronic  and  isobaric  analogue,  F  .  Thus, 
in  absence  of  data  on  NaOH,  we  can  use  available  information 
on  NaF  and  other  alkali  halides  to  try  to  understand  the  properties 
of  the  alkali  hydroxides.  There  have  been  a  number  of  molecular 
beam  scattering  measurements  of  alkali  halides  with  other  alkali 
halides  or  hydrogen  halide  molecules.43  As  a  class,  they  all  imply 
the  existence  of  a  strong  collision  complex.  In  the  case  of  CsCl 
+  KI  -*  Csl  +  KC1,  the  reaction  proceeds  without  any  activation 
energy  (despite  being  four-centered)  and  has  a  large  total  cross 
section.4*  This  bas  been  explained  by  the  large  ionic  character 
of  both  species  and  by  the  fact  that  alkali  halide  dimers  have  large 
binding  energies  in  the  geometry  of  a  cyclic  planar  rhomboid.43 
This  ion-pair  intermediate  formulation  has  also  been  invoked  in 
describing  the  Bal}  +  HCI  reaction,41  where,  for  large  alkaline 
earths  such  as  barium,  we  can  treat  Bal2  as  (Bal)+I".  What  we 
can  draw  from  these  comparative  systems  is  that  a  gas  kinetic 
rate  constant  without  any  activation  barrier  is  a  reasonable  ex¬ 
pectation  for  this  system.  In  the  extrapolation  of  the  measured 
rate  data  to  lower  temperatures  characteristic  of  the  stratosphere, 
one  would  not  expect  any  great  dependence  on  temperature  for 
the  reaction  rate. 

The  importance  of  reaction  2  in  stratospheric  ozone  chemistry 
may  be  estimated  by  comparing  the  Cl  regeneration  rates  from 
HCI  via  reaction  with  OH  with  regeneration  rates  from  reaction 
with  NaOH  followed  by  photolysis  of  NaCl.  Estimates  of 
stratospheric  NaCl  photolysis  rates  by  Rowland  and  Rogers24  are 
in  the  range  from  10"2  to  10"3  s"1.  At  40  km,  the  meteoric  sodium 
is  partitioned  among  the  species  NaOH,  NaO,,  and  NaO.  Es¬ 
timates  of  total  stratospheric  sodium  concentrations  by  Liu  and 
Reed10  are  on  the  order  of  5  X  103  cm"3.  If  all  of  the  meteoric 
sodium  were  in  the  form  of  NaOH  at  40  km,  our  measurement 
of  k 2  would  give  a  first-order  rate  constant  for  Cl  formation  of 
1.4  X  10~*  s"',  providing  that  /N<C1  is  in  the  range  estimated  by 
Rowland  and  Rogers.  This  is  a  factor  of  —20  faster  than  the 
first-order  rate  constant  for  Cl  regeneration  from  OH  +  HCI  of 
6  X  10-4  s'1,  based  on  an  OH  concentration  of  107  cm"3  and  the 
value  k5(250  K)  ■  5.7  X  10"' 3  cm3  molecule"1  s'*  .*'  Although 
the  partitioning  of  total  sodium  among  NaOH,  NaO,  and  NaO, 
requires  a  more  detailed  knowledge  of  rate  constants  for  the 
processes  shown  in  Figure  1,  this  calculation  demonstrates  that, 
even  if  only  5%  is  in  the  form  of  NaOH,  regeneration  of  Cl  from 
HCI  via  alkali  chemistry  would  be  comparable  to  regeneration 
by  the  OH  reaction. 

This  simple  computation  leads  us  to  two  conclusions  First, 
our  results  for  this  and  the  other  alkali  reactions  described  show 
that  meteoric  metal  reactions  may  have  a  potentially  significant 
impact  on  our  understanding  of  chemistry  in  the  mesosphere  and 
upper  stratosphere.  The  large  values  of  the  sodium  rate  constants 
measured  in  this  study  emphasize  this  possibility  Secondly,  it 
clearly  motivates  the  need  for  further  investigation  of  this  chem- 
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istry,  via  detailed  models,  additional  laboratory  measurements  of 
reaction  rate  constants  and  photolysis  cross  sections,  and  direct 
measurement  of  molecular  sodium  concentration  profiles  in  the 
upper  atmosphere. 
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The  gaseous  diffusion  coefficients  for  sodium  and  potassium  atoms  have  been  measured  as  a 
function  of  temperature  in  a  flow  tube.  Analysis  is  based  on  a  simple  relationship  between 
diffusion  and  the  observed  alkali  decay  rate  for  the  case  of  very  reactive  walls.  The  resulting  values 


for  the  diffusion  coefficient  in  units  of  cm2  s 

1  at  1  atm  are: 

N*-He 

(2.2  ±  0.3)X  10"’ 

309-473  * 

-■ 

Na-N2 

(8.o±o.g)xio -*r”,iain 

320-698  K 

• 

Ne-Ar 

|1.2±0.9)XKr57",,ta‘' 

322-350  K 

K-He 

(1.3  ±  1.01x10—  T'u±0" 

301-734  K 

k-n2 

(4.5  ±  2.8)  x  10-5  7” 541009 

302-720  K 

The  results  agree  well  with  previous  measurements  and  the  comparison  with  theoretical 
expectations  is  discussed. 


I.  INTRODUCTION 

The  transport  properties  of  gases  are  important  in  in¬ 
vestigating  the  details  of  atomic  and  molecular  interaction 
potentials  and  in  a  more  practical  sense,  are  needed  to  de¬ 
scribe  the  dynamics  of  flow  and  reaction  in  natural  and  in¬ 
dustrial  processes.  Our  interest  in  diffusion  coefficients 
stems  from  the  need  to  correct  flow  tube  kinetic  rate  data  for 
wall  and  diffusional  effects. 

Measurements  of  diffusion  coefficients  of  alkali  atoms 
date  back  to  the  early  1 930’s.  Since  then,  a  variety  of  experi¬ 
mental  approaches  have  been  used,  including  static  diffusion 
cells,'  spin  relaxation,2'5  pulsed  irradiation,5-®  flame  pho¬ 
tometry,9-12  resonance  ionization  spectroscopy,1 3-14  and  so¬ 
phisticated  digital  correlation  techniques.15  The  results  re¬ 
ported  in  this  paper  are  obtained  from  the  observed  decay  of 
laser  induced  alkali  atom  fluorescence  in  a  flowing  gas  in  a 
regime  where  wall  loss  is  very  efficient  and  axial  diffusion 
unimportant.  The  techniques  discussed  in  analyzing  the  data 
are  a  special  case  of  the  more  general  solution  used  in  cor¬ 
recting  observed  bimolecular  reaction  rate  constants  for  the 
effects  of  wall  removal  and  diffusion  in  a  flow  tube. 16-19  In 
this  work,  the  bimolecular  reaction  rate  is  zero,  and  the  solu¬ 
tion  is  somewhat  simplified. 

II.  EXPERIMENTAL 
A.  Apparatus 

The  high  temperature  fast  flow  reactor  used  in  these 
experiments  has  been  described  in  detail  previously.20  Brief¬ 
ly,  a  7.26  cm  diameter,  120  cm  long  alumina  tube  is  heated 
with  Kanthal  resistance  heaters  capable  of  maintaining  uni¬ 
form  temperatures  in  the  reaction  zone.  The  helium,  argon 
or  nitrogen  carrier  gas  enters  the  flow  tube  through  two  mul- 
lite  multichannel  arrays  which  both  preheat  and  laminarize 
the  flow  upstream  of  the  reaction  zone.  A  sufficient  distance 
is  allowed  downstream  of  these  arrays  for  the  carrier  gas  to 
develop  a  parabolic  flow  profile  before  reaching  the  reaction 
zone.  The  gas  is  pumped  by  a  450  1  s'  1  Roots  blower  and 
mechanical  fore  pump  combination.  Linear  flow  velocities 
from  220  to  2800  cm  s  1  are  used  in  the  present  experiments. 


Carrier  gas  flow  rates  from  30  to  400  STP  cm3  s' 1  are  mea¬ 
sured  with  calibrated  rotameters.  Total  pressure  is  measured 
with  a  capacitance  manometer  at  the  downstream  end  of  the 
reaction  zone.  Gas  temperatures  are  obtained  with  a  chro- 
mel-alumel  thermocouple  which  is  movable  throughout  the 
reaction  zone.  Although  the  apparatus  is  capable  of  tern 
peratures  up  to  1500  K,  experimental  considerations  as  de 
scribed  below  limited  the  upper  temperature  to  700  K  for 
these  studies. 

For  nominal  room  temperature  measurements  (flow 
tube  wall  heaters  off),  the  axial  temperature  profile  is  noi 
uniform,  but  exhibits  a  small  drop  in  temperature  down 
stream  of  the  alkali  oven.  Gas  passing  by  or  through  the  oven 
is  heated,  and  then  cools  as  it  travels  further  down  the  tube. 
Over  the  region  in  which  decay  measurements  are  taken,  the 
maximum  variation  from  the  average  temperature  is  typical¬ 
ly  five  to  ten  degrees,  depending  on  the  flow  velocity,  total 
pressure,  and  identity  of  the  carrier  gas.  The  more  importam 
effect  is  that  the  average  gas  temperature  is  1 5’  to  25*  warmer 
than  room  temperature,  and  is  the  correct  value  to  use  in 
analyzing  the  data.  At  elevated  temperatures  ( >  400  K),  gas 
heating  is  dominated  by  the  walls  and  mullite  arrays,  and  the 
oven  contribution  is  secondary.  As  a  result,  the  axial  profiles 
are  more  uniform,  with  variations  of  only  2*-5*. 

Alkali  atoms  are  generated  by  heating  the  pure  metal  in 
a  25  mm  diam  cylindrical  silver  plated  monel  oven  to  a  tem¬ 
perature  sufficient  to  obtain  a  vapor  pressure  of  1 0  “  6  to  1 0  ~ 4 
Torr  within  the  oven.  The  oven  is  mounted  on  the  end  of  a 
movable  13  mm  o.d.  alumina  tube  concentric  to  the  main 
flow  tube.  The  vapor  is  entrained  in  a  flow  of  carrier  gas  and 
introduced  into  the  flow  tube  in  one  of  two  configurations  In 
the  first,  the  oven  is  placed  at  a  fixed  position  upstream  of  the 
multichannel  arrays.  The  alkali  vapor  is  introduced  into  the 
main  flow  through  a  10  cm  length  of  19  mm  o.d.  silver  tubing 
which  passes  through  the  center  of  the  arrays.  In  the  second 
configuration,  the  oven  is  placed  downstream  of  the  arrays 
and  is  movable  throughout  the  reaction  zone.  The  oven  tem¬ 
perature  is  controlled  with  resistive  heating  elements  inde¬ 
pendent  of  the  flow  tube  heaters.  Initial  alkali  atom  concen¬ 
trations  in  the  reaction  zone  are  maintained  at  less  than  10"’ 


J.  Chem  Phys.  81  (11).  1  December  1984  0021  -9606/84/2351 25-06S02  10 


.  c  1984  American  Institute  ot  Physics 


61 


5t26 


Joel  A.  Silver:  Atomic  Na  and  K  diffusion  coefficients 


Residence  Time  (sec) 

FIG  1 .  T  ypicaJ  decay  plots  of  Na  using  fixed  and  moveable  source  methods. 
Both  runs  were  made  at  469  K  in  1  02  Torr  N,  |al  Movable  source; 
(u>  =  1080cms_l,2=  1 5-65  cm  (b|  Fixed  source:  z  »  65cm,  (a)  =  850- 
2750 cm  s~ 

cm  ”  *  by  adjusting  the  oven  temperature  and  the  carrier  gas 
flow  rate  through  the  oven. 

The  alkali  atoms  are  detected  by  laser  induced  fluores¬ 
cence  at  the  downstream  end  of  the  flow  tube.  A  Molectron 
DL14  nitrogen  pumped  dye  laser  is  used  to  excite  the  4s 
251/2— *5p  2PU2  transition  at  404.4  nm  for  potassium  and  the 
is  2Sin—*ip  2P,/2  transition  at  589.0  nm  for  sodium.  Flu¬ 
orescence  is  collected  using  a  gated  phototube  and  a  comput¬ 
er  controlled  data  acquisition  system. 20  The  signals  are  aver¬ 
aged  over  100  laser  pulses  after  subtracting  nonfluorescent 
background  contributions  and  normalizing  for  pulse  to 
pulse  fluctuations  in  laser  intensity.  To  verify  that  the  transi¬ 
tions  were  not  saturated  and  that  the  samples  were  optically 
thin,  it  was  observed  that  the  measured  fluorescence  varied 
linearly  with  both  alkali  number  density  and  laser  intensity. 
Although  no  direct  calibration  of  the  fluorescence  signal  was 
attempted,  estimates  of  sensitivity  from  known  phototube 
and  integrator  responses,  measured  laser  power,  and  atomic 
transition  probabilities  indicate  a  detection  limn  for  Na  of 
104  cm  ”  ’  and  for  K.  of  106  cm  1  using  these  transitions. 

Loss  of  atoms  on  the  walls  as  a  function  of  residence 
time  in  the  flow  tube  is  measured  in  two  ways.  The  oven 
position  z  is  varied  at  fixed  average  flow  velocity  ( u ).  The 
observed  decay  rate  k  *  is  equal  to  the  slope  of  In  (signal]  vs 
time(z/<u)i  The  second  method  uses  the  fixed  alkali  source. 
Here  the  residence  time  .s  varied  by  changing  the  flow  veloc¬ 
ity,  with  the  pressure  constant,  so  that  the  diffusion  coeffi¬ 
cient  is  the  same  for  all  measurements  A  plot  of  In  (sig¬ 
nal-  (u)\  vs  time  gives  a  slope  equal  to  k  •  The  velocity  term 
in  the  ordinate  is  required  to  compensate  for  the  change  in 
alkali  density  in  the  detector  due  to  dilution  by  the  increased 


carrier  gas  flow  as  the  velocity  is  changed.  This  latter  tech¬ 
nique  is  often  useful  when  it  is  difficult  to  use  a  movable 
source.  Comparison  of  the  two  methods  shows  no  significant 
difference  in  the  final  observed  rate  constants  under  the  con¬ 
ditions  present  in  this  work.  An  example  of  results  for  both 
types  of  measurements  under  the  same  pressure  and  tem¬ 
perature  conditions  is  shown  in  Fig.  1 .  Decays  are  very  linear 
and  exhibit  little  scatter.  Most  of  the  measurements  for  Na 
were  made  using  the  variable  distance  method,  while  all  of 
the  K  runs  used  a  fixed  distance  of  73  cm.  At  temperatures 
above  700  K,  alkali  atoms  were  observed  to  effuse  from  the 
walls,  as  evidenced  by  nonzero  signals  with  the  oven  entrain¬ 
ment  flow  turned  off.  As  a  result,  no  data  were  taken  above 
this  temperature. 


B.  Analysis 


Data  analysis  is  based  on  the  solution  to  the  diffusion 
equation  for  laminar,  Poiseuille  flow  in  a  cylindrical  tube  in 
which  the  only  loss  process  is  by  heterogeneous  removal  at 
the  wall. I<u8  This  equation  is  expressed  in  terms  of  the  tube 
radius  a,  radial  distance  r,  axial  distance  z,  species  concen¬ 
tration  C,  diffusion  coefficient  De,  and  mean  velocity  (u)  as 


,  r*\dC  _  (PC  1  <9C  PC\ 


m 


Note  that  this  contains  both  axial  and  radial  diffusion  terms. 
The  first  order  rate  constant  for  wall  removal  kw  appears  in 
the  boundary  condition 


(2) 


These  equations  can  be  rewritten  in  terms  of  dimensionless 
parameters  as 


(1  R  ]lz~D \1P+  RdK  +  dz2 T 


(3) 


with 


(4) 


where  R  =  r/a,  Z  -  z/a,  D  —  Dc/2a{u),  and  Kw  —  kwa/ 
2{u).  Furthermore,  y,  the  collision  efficiency  for  wall  loss  of 
C  (also  called  the  sticking  coefficient),  is  related  to  kw  by 
elementary  kinetic  theory21 


k  SL 

w  2a(\-y/2)' 


15) 


where  c  is  the  mean  molecular  velocity  of  C. 

Although  the  exact  solution  of  Eq.  (3),  subjected  to  the 
boundary  conditions,  can  be  obtained  only  by  numeric 
methods,  it  can  be  solved16”18  by  assuming  a  solution  of  the 
form 


C=  £  A,g,{R)e-KTZ  .  (61 

Substitution  of  this  into  Eqs.  (3)  and  (4)  leads  to  an  infinite  set 
of  ordinary  differential  equations 

dzg,(R  )  +  1  dg,(R  1 
dR 2  R  dR 

+  [K?1  +  (\-R2)KyD]gl[R)  =  0  (7) 
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FIG.  2.  Observed  decay  rate  [k  •)  vs  sticking 
coefficient  for  a  series  of  diffusion  coefficients. 
The  dashed  line  corresponds  to  the  situation 
where  axiaJ  and  radial  effects  are  equal.  The 
symbols  (x)  on  the  right-hand  axis  corre¬ 
sponds  to  k*  calculated  using  Eq  (15). 


Sticking  Coefficient  (y) 


<&(*>]  _  _  *u ,*■(*!*., 

.  dR  J*.,  2D  (  ’ 

It  is  found  that  g,  (R  )  can  be  expressed  as  a  series  in  even 
powers  of  R  (dropping  the  subscripts  for  clarity) 


dR  )  =  £  BnR  *  , 

n  —  0 

with  coefficients 


fl,= 

and 


aBq 
4  ’ 


where 


(9) 

(10) 

(U) 

(12) 


a  =  **J  +  lT  (13> 

Substitution  of  Eq.  (9)  into  the  boundary  condition  (8)  leads 
to  the  expression  (for  each  0 

£  Bn  (ln  +  ^  =  F[Kwj:*J))=Q,  (14) 

n  -  o  \  UJ  / 

If  K,m  and  Dare  specified,  then  the  first  positive  root  oTF 
is  K  the  second  root  is  K  *,  and  higher  roots  become  K 
K  J,  etc.  For  sufficiently  large  values  of  2,  the  higher  order 
terms  in  Eq.  (6)  rapidly  decay,  so  that  the  observed  disap¬ 
pearance  of  C  is  described  by  a  single  exponential  term 
W  =  1),  with  the  corresponding  decay  rate  K  *  related  to  the 
observed  decay  rate  k*  (s'1)  such  that  AT*  =  ak  J/(u). 

Figure  2  illustrates  the  relationship  between  K  *  (ob¬ 
served  wall  loss  rate  constant!  and  the  sticking  coefficient  y 


as  a  function  of  the  diffusion  parameter  D  for  sodium  at  300 
K.  As  shown  in  this  figure,  for  smaller  values  of  D,  the  K  * 
curve  flattens  as  y  approaches  unity.  For  y>0. 1,  K*  is  inde 
pendent  of  y  and  varies  linearly  with  D.  These  phenomena 


1/P  (torr  ’) 


FIG.  3  First-order  decay  rate  Ik;  I  vs  \/P  for  Na  in  N.  at  three  tempera 
lures. 
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Temperature  (K) 

300  «00  600  MO  1400  2000  2600 


FIG  4  Diffusion  coefficients  for  sodium  as  a  function  of  temperature.  Ex¬ 
perimental  data  in  helium  are  shown  for  this  work  (O),  Fairbank  el  al.  (Ref. 
151 1 VI.  Husain  and  Plane  (Ref.  6)  (•),  Ramsey  and  Anderson  (Ref.  2)  (CH), 
and  Hanel  era/  (Ref.  !)(<)),  and  Bichi  eiai  (Ref.  3)(0);  in  nitrogen  for  this 
work  I  Aj.  Husain  and  Plane  (Ref.  6)  (A).  Ramsey  and  Anderson  (Ref.  2)  (■), 
and  Handera/  (Ref  1 1  (41,  m  argon  for  this  work  |  ☆),  Fairbank  rial.  I  Ref. 
151  (T|,  and  Hand  er  al.  (Ref  1)  (*);  in  flames  for  Ashton  and  Hayhurst 
iRef.  10)  (A)  and  Snelleman  (Ref.  12)(  +  |.The  crosses  (x)  correspond  to 
calculations  of  Redko  (Ref.  26)  in  helium,  and  the  solid  lines  the  best  fits  to 
the  present  data.  The  dashed  line  is  an  extrapolation  of  the  fit  in  N,  to  flame 
temperatures. 

occur  when  the  decay  of  C  is  controlled  by  ra dial  diffusion. 
The  effects  of  axial  and  radial  diffusion  become  equal  at 
D  —  0.072  and  axial  diffusion  dominates  at  larger  values  of 
D.  Ferguson  et  al. 19  have  shown  that  in  the  limit  of  unit  wall 
loss  [y=  I)  and  radial  diffusion  control, 

D  =  0.137  K*.  (15) 

The  symbols  (X )  on  the  right-hand  axis  of  Fig.  2  represent 
K  *  calculated  for  the  given  values  of  D  using  this  formula. 
Agreement  between  these  and  the  exact  solutions  are  quite 
good  in  the  radial  diffusion  limit  (Z><0.05). 


III.  RESULTS 

A  total  of  29  wall  loss  measurements  were  made  for 
sodium  and  eight  for  potassium  between  300  and  700  K.  At 
each  temperature,  runs  were  made  at  a  number  of  different 
pressures  in  the  range  of  1  -8  Torr.  All  runs  were  in  the  radial 
diffusion  limit.  Since  alkali  atoms  are  removed  by  the  un¬ 
coated  alumina  flow  tube  walls  with  near  unit  efficiency, 
diffusion  coefficients  were  directly  obtained  from  the  ob¬ 
served  decay  rates  using  a  computer  program  based  on  the 
derivation  in  Sec.  II  B  It  is  assumed  that  y  —  I,  but  even  if 
y  —  0.2,  the  error  in  D,  would  be  <  5%.  The  assumption  of 
y  =  1  is  further  confirmed  by  observing  that  since  Dc  is  in¬ 
versely  proportional  to  pressure,  a  plot  of  k  J  vs  \/P  for  a 
given  buffer  gas  and  temperature  will  be  linear  only  in  the 
limit  of  y— *  I  An  example  of  this  linearity  is  seen  in  Fig,  3  for 
Na  in  N2  at  three  different  temperatures. 

As  mentioned  in  Sec.  II  B,  it  was  assumed  that  Poi- 


T tmptrtlur*  (K) 


300  400  600  900  1400  2000  2600 


FIG  5.  Diffusion  coefficients  for  potassium  as  a  function  of  temperature. 
Experimental  data  in  helium  are  shown  for  this  work  (O),  Bemheun  and 
Korte  (Ref.  5)  |+),  and  Husain  and  Plane  (Ref.  7)  {•);  in  nitrogen  for  this 
work  (/\)  and  Husain  and  Plane  (Ref.  7)  (A);  in  flames  for  Hayhurst  and 
Ashton  (Ref.  1 1)  (*).  The  crosses  I X I  correspond  to  calculations  of  Redko 
(Ref.  26)  in  helium  and  the  solid  lines  the  best  fits  to  the  present  data.  The 
dashed  line  is  an  extrapolation  of  the  fit  in  N,  to  flame  temperatures. 


scuffle  flow  exists  in  the  measurement  region  and  that  only 
the  first  term  of  Eq.  (6)  is  significant.  For  parabolic  flow  to  be 
fully  developed,  an  inlet  mixing  length  of  approximately 
R,d/ 20  is  required,22  where  Rr  is  the  flow  Reynold’s  num¬ 
ber  and  d  the  tube  diameter.  In  ail  experiments  this  distance 
is  <  20  cm,  permitting  a  minimum  of  65  cm  for  measuring 
decay  rates. 

The  second  requirement  is  that  the  decay  of  C  be  de¬ 
scribed  by  a  single  exponential  term.  This  requires  calcula¬ 
tion  of  the  A,  andg,(R  )  terms  in  Eq.  (6).  Following  the  proce¬ 
dure  of  Pirkle  and  Sigillito, 17  the  first  eight  terms  of  C{rj) 
were  calculated.  Initial  conditions  accurately  reflected  the 
true  source  inlet  geometry,  rather  than  assuming  an  initial 
uniform  concentration  profile  of  C  across  the  entire  flow 
tube  as  done  by  Pirkle.  The  results  show  that  the  sum  of  all 
higher  order  modes  decay  to  less  than  5%  of  the  first  term  in 
less  than  10  cm.  Furthermore,  no  curvature  was  observed  in 
the  decay  data  at  shorter  reaction  times  as  might  be  expected 
if  higher  order  terms  were  significant  (Fig.  1). 

Diffusion  coefficients  resulting  from  these  experiments 
are  shown  in  Fig.  4  for  sodium  and  Fig.  5  for  potassium.  The 
data  (including  their  associated  uncertainties!  have  been  fit 
to  a  functional  form  Dc  =  aT"  using  a  nonlinear  least  square 
analysis.  The  best  fit  coefficients  are  shown  in  Table  I  and  the 
solid  lines  in  Figs.  4  and  5  represent  these  fits.  The  expressed 
total  experimental  uncertainties  for  each  Dc,  including 
allowance  for  systematic  errors,  can  be  estimated  as  the 
square  root  of  the  sum  of  the  squared  individual  uncertain¬ 
ties  due  to  (a|  flow  velocity,  temp'rature.  and  pressure  cali¬ 
bration  factors,  7.2%,  (b)  the  precision  in  determining  the 
coefficients  a  and  n  in  the  expression  for  Dc.  5.4%,  and  (cl 
the  uncertainty  in  converting  K  *  to  Dc,  5%. 
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TABLE  I.  Comparison  of  experimental  values  of  diffusion  coefficients  for  sodium  and  potassium. 


Alkali 

Buffer  gas 

m> 

D,  (cm2  s" 1  at  1  atm;  ±  la) 

Method 

Reference 

Na 

He 

304-473 

(2.2  ±0.3)x!0-57"  ” 1002 

Flow  tube 

Thu  work 

N, 

320-698 

(8.0±0.8|xl0-‘ 

Ar 

322-350 

(1.2  ±  0.9)X  10~5  rl,lt0" 

Na 

He 

724,844 

3.58  ±  0.09,  1.36  ±0.11  (2*) 

Pulse  irradiation/ 

Husain  and 

N, 

724,844 

1.17  ±0.04,  1.14  ±0.05  (2 a) 

resonance  absorption 

Plane  I  Ref.  6| 

CO, 

724,844 

4.96  ±  0.04,  1.63  ±0.05  (2*) 

Na 

He 

300 

0.474  ±  0. 120,  0.564  ±  0.126  (3o) 

Laser  beam  transit 

Fairbank 

Ar 

300 

0.238  ±  0.053, 0.329  ±  0.080  (3<r) 

times 

etal.  (Ref.  15) 

Na 

He 

453 

1.2  ±0.2 

Spin  relaxation 

Bicchi  ei  al  (Ref.  3| 

Ne 

453 

0.67  ±  0.05 

!°Na 

Ne 

300 

0  35  ±  0.08 

Pulsed  irradiation/ 

Coolen  and 

- 

resonance  fluorescence 

Hagedoom  (Ref  8) 

Na 

H, 

473 

1.3 

Spin  relaxation 

Ramsey  and 

He 

428 

0.92 

Anderson  (Ref.  2) 

Ne 

428 

0.50 

N, 

453 

0.54 

Na 

H, 

655 

3.14 

Static  diffusion 

Hartel  et  al.  (Ref  1) 

He 

654 

2.17 

cel) 

N, 

6JJ 

0.91 

Ar 

654 

0.88 

QH„ 

655 

0.23 

Na 

C,H,/0;/N, 

Flame  products 

2440 

9.9  ±  0.5 

Flame  photometry- 

Snelleman  (Ref  12) 

Na 

H,/0,/N, 

1920-2520 

9.3  x  io- 7  r11 

Flame  photometry" 

Ashton  and 

Flame  products 

Hayhurst  IRef  101 

K 

He 

301-724 

(i.s±i.o)xio-4r,a*0'’ 

Flow  lube 

Thu  work 

N, 

302-720 

(4.5±2.8tXlO-’r',‘t00’ 

K. 

He 

753,873 

6.41  ±0.41, 4.97  ±0.34(2*) 

Pulse  irradiation/ 

Husain  and 

N, 

753,873 

2.70  ±0.18, 1.54  ±0.11(2*) 

resonance  absorption 

Plane  (Ref.  7) 

CO, 

753,873 

4.49  ±0.41,  1.37  ±0.34(2<x> 

K 

H,/0,/N, 

1900-2520 

3.ixio-7r“ 

Flame  photometry' 

Ashton  and 

Flame  products 

Hayhurst  (Ref.  11) 

K 

He 

358 

1.2 

Spin  relaxation 

Bernbeim  and 

Korte  (Ref.  5) 

‘  Ashton  reports  the  results  of  Snelleman  to  have  if”  temperature  dependence  although  no  range  of  temperatures  is  specified. 
"Average  of  eight  data  points  At  to  an  expression  of  the  form  aT" .  Individual  data  are  shown  m  Fig.  4. 

*  Average  of  eight  data  points  fit  to  an  expression  of  the  form  aT" .  Individual  data  are  shown  in  Fig.  3. 


IV.  DISCUSSION 

The  results  presented  here  compare  well  with  previous 
measurements  for  sodium  and  potassium  (Table  I)  and  span 
a  wide  temperature  range,  permitting  an  accurate  determin¬ 
ation  of  the  temperature  dependence  of  Dc .  Despite  the  var¬ 
iety  of  different  techniques  used,  the  data  agree  well  over  a 
2000  K  temperature  range,  based  on  an  extrapolation  of  our 
temperature  dependences.  More  rece  /:»asurements  in¬ 
clude  those  of  Ashton  and  Hayhurst, '  <vho  measured  the 
diffusion  of  alkali  atoms  (Na,  K,  Cs,  Rb)  and  LiOH  in  H:/ 
O./Nj  flames,  introduced  from  a  point  source.  The  data  for 
Na  and  K  match  the  extrapolated  N2  fits  very  well,  as  might 
be  expected  since  the  flame  products  are  54%-75%  N2,  the 
balance  being  mostly  H20  with  some  (1%-10%|  Hj*  Hu¬ 


sain  and  Plane6,7  obtained  diffusion  coefficients  in  a  pulsed 
irradiation  cell,  reporting  substantially  different  results  at 
the  two  temperatures  used.  Their  lower  temperature  data 
(724  and  753  K  for  Na  and  K,  respectively)  appear  to  be 
anomalously  high,  but  no  reason  for  this  is  given.  Fairbank  et 
al. , 15  in  contrast  to  this  and  all  of  the  other  work,  use  a  much 
different  approach  to  measure  Dc  Transit  time  measure¬ 
ments  of  a  single  atom  across  a  laser  beam  are  computed 
from  the  observed  photon  burst  signal.  In  these  experiments, 
the  atom  absorbs  and  reemits  the  resonance  laser  radiation 
many  times  and  the  time  dependent  signal  is  tre  .  J  hy  digi¬ 
tal  correlation  techniques  to  provide  the  mean  traversal  time 
and  thus  the  diffusion  coefficient. 

The  calculation  of  Dc  from  transport  theory  has  been 
well  described.21  Given  the  interaction  potential  energy  sur- 
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face  for  the  diffusing  species  and  the  buffer  gas,  Dc  can  readi¬ 
ly  be  calculated.  In  the  absence  of  experimental  data,  ap¬ 
proximations  using  simplified  potentials  such  as  a 
Lennard- Jones  often  lead  to  values  of  Dc  accurate  to  within  a 
factor  of  2.  For  the  alkali  atom-rare  gas  systems,  more  de¬ 
tailed  interaction  potentials  are  now  available, :<  M  and 
Redko26-27  has  used  them  to  calculate  diffusion  integrals  and 
coefficients  for  ground  and  excited  electronic  states  of  Na 
and  K  in  rare  gases  over  the  temperature  range  300  to  3000 
K.  His  results  for  Na(3  2S )  and  K(4  2S )  in  He  are  shown  as 
crosses  (jc)  in  Figs.  4  and  5,  respectively,  and  show  remarka¬ 
bly  good  agreement  with  the  experimental  results  over  the 
temperature  range  studied,  including  the  different  tempera¬ 
ture  dependence  of  Dc  between  Na  and  K. 

The  measurements  of  diffusion  coefficients  in  our  labo¬ 
ratory  are  presently  being  extended  to  include  molecular  al¬ 
kali  species.  Preliminary  results  for  NaOH  in  He  at  308  K 
are  0.47  ±  0.05  and  0.4  ±  0. 1  cm2  s" 1  for  NaCl  in  He  at  1 
atm. 
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